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I.  INTRODUCTION 

The  Advanced  Sensors  Directorate,  U.S.  Army  Missile  Research  and 
Development  Command  (MIRADCOM),  has  conducted  an  extensive  exploratory 
development  program  In  Imaging  seekers  for  application  on  terminal 
homing  missiles.  Several  prototype  seekers  have  been  produced  featuring 
a variety  of  sensors,  detector  arrays,  and  scanning  mechanisms.  However, 
the  ability  to  characterize  the  effects  of  the  atmosphere  on  this  class 
of  seekers  Is  not  well  established,  consequently,  predicting  performance 
of  a terminal  homing  weapon  using  this  seeker  under  a variety  of  meterol- 
oglcal  conditions  has  not  been  possible. 

The  general  objective  for  this  effort  was  to  develop  the  analytical 
description  of  the  atmospheric  effects  on  an  imaging  seeker  that  allows 
seeker  performance  to  be  related  to  the  specific  meteorological  conditions 
over  the  line  of  sight  path  from  the  seeker  to  the  target.  Inclusion  of 
detailed  atmospheric  effects  in  the  parametric  analysis  of  sensor  perfor- 
mance produces  added  complexity  which  is  not  in  every  case  separable  from 
effects  such  as  the  modulation  transfer  function  (MTF)  of  the  imaging 
optics.  Further,  the  theoretical  development  for  absorption  and  scattering 
effects  is  provided  quite  distinctly  from  that  for  the  effects  of 
refractive  Index  fluctuation.  Contrast  reduction  due  to  scattering  and 
absorption  is  based  on  thermodynamic  arguments  and,  if  certain  assumptions 
are  accepted.  Is  derived  In  a form  which  Is  readily  compared  with 

measured  meteorological  parameters.  On  the  other  hand,  the  effects  of 

» 

atmospheric  turbulence,  or  refractive  Index  fluctuation,  is  based  on  one 
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of  several  approximate  solutions  to  Maxwell's  equations  and  results  in 
an  expression  dependent  on  an  ensemble-averaged  "wave  structure  function" 
for  the  atmospheric  path. 


It  was  desired. to  develop  an  analytical  description  of  the  effects 
of  the  atmosphere  on  imaging  seekers  which  includes  scattering,  absorption 
and  refraction  effects  in  a single  theoretical  development.  Specifically, 
the  theory  should  show  the  effects  of  target  size  on  contrast  transfer 
through  the  atmosphere,  include  spectral  dependence  with  an  analysis  of 
the  expected  effects  of  finite  spectral  band  imaging  and  provide  the 
basis  for  the  development  of  an  atmospheric  measurements  program. 

The  major  body  of  knowledge  about  the  effects  of  the  atmosphere 
on  optical  radiation  is  best  described  by  a small  number  of  primary 
references.  One  goal  of  the  current  effort  was  to  establish  a list  of 
these  primary  references  and  to  present  in  easily  readable  summary  form, 
the  specific  atmospheric  effects  described.  An  alternative  objective  was 
to  identify  deficiencies  in  the  theoretical  development  as  related 
specifically  to  imaging  seeker  applications.  This  was  accomplished 
through  comparison  of  assumptions  and  approximations  made  in  the  theoreti- 
cal development  with  expected  conditions  for  the  seeker. 


Specific  casks  Investigated  and  the  results  are  Identified 


[below. 

Task  a)  Develop  theory  for  the  purpose  of  treating  the 
atmosphere  between  an  Imaging  device  and  an  ob- 
ject as  an  optical  window  to  which  techniques 
for  quantitative  description  such  as  the  modul- 
ation transfer  function  can  be  applied.  Struc- 
ture theory  Into  a mathematical  model  defining 
key  parameters  which  must  be  known. 

Results  of  Task  a) 

Starting  with  Maxwell's  equations  (32) -(35)  and  proceeding  through 
a series  of  assumptions  and  approximations  we  arrive  at  equation  (73), 
the  scalar  Helmholtz  wave  equation,  which  describes  optical  propagation 
in  the  atmosphere.  Several  solution  methods  to  this  second-order 
differential  equation  are  examined.  A rigorous  treatment  based  on 
Green's  theorem  is  described  in  Section  3.2.4,  yielding  equation  (130), 
which  describes  the  optical  amplitude  function  at  i • d in  terms  of  the 
input  optical  amplitude  at  z » 0 and  atmosphere  dependent  parameters. 

Section  4.1  extends  the  result  to  give  intensity  in  the  image 
plane  of  a simple  imaging  system.  From  this  is  computed  the  atmosphere- 
lens  point  spread  tunctlon  (equation  157),  the  optical  transfer  function 
(equation  168)  and  the  MTF  (equation  169). 

Task  b)  Develop  theory  for  comparing  atmospheric  effects  on 


the  modulation  transfer  function  (MTF)  and  on  the 
reduction  of  contrast.  Specifically  examine  the 


effects  of  spatial  frequency  on  the  reduction 
of  contrast  by  the  atmosphere  and  relate  the 
theories  for  MTF  and  contrast  transfer  function 
(CTF) . 

Results  of  Task  b) 

Contrast  based  on  the  model  developed  in  task  a)  is  computed  which 
includes  dependence  on  spatial  frequency  content  of  the  target  object, 
atmospheric  parameters  and  ambient  irradiance  in  Section  4.4.  Universal 
contrast  and  modulation  contrast  through  the  atmosphere  are  given  for 
a periodic  bar  chart  in  equations  (190)  and  (192).  Appendix  C gives 
further  details  on  MTF  and  CTF. 

Task  c)  Develop  a mathematical  model  for  treating  the 
effects  of  the  atmosphere  for  the  finite  spec- 
tral band  application  (i.e.  non-monochromatic) . 

Emphasis  will  be  on  the  visible  and  near  infra- 
red spectral  regions  (0.4  to  1.2  micrometers). 

Results  of  Task  c) 

The  theory  of  Duntley  and  Middleton  for  contrast  reduction  is 
extended  to  the  finite  spectral  band  case  in  Section  2.2  by  defining 
a broadband  extinction  coefficient.  Atmospheric  transmittance  is  then 
compared  for  the  new  theory  versus  transmittance  at  a central  wavelength. 
The  effor  is  small  except  for  marginal  conditions  (i.e.,  slant  range 
near  the  visibility  range). 


Task  d)  Develop,  where  needed,  techniques  for  measuring 


the  parameters  required  In  an  atmospheric  MTF 
model.  Recommend  hardware  which  presently  exists 
or  hardware  which  must  be  developed  for  measuring 
the  desired  parameters. 

Results  of  Task  d) 

Parameter  measurement  methods  are  discussed  in  Section  5.0  for 
atmospheric  extinction  coefficient  and  atmospheric  refractive  index 
structure  function.  Additionally,  measurement  methods  for  determining 
spatial  size  dependence  on  contrast  transfer  are  discussed  briefly. 

This  effort  did  not  get  to  the  point  of  defining  measurement  hardware. 
Additional  measurement  program  planning  is  needed. 

Task  e)  Having  developed  and  proven  validity  of  an 

atmospheric  MTF  model  in  the  visual  region  of 
the  light  spectrum,  extend  model  to  encompass 
the  near  infrared  (0.7  to  0.9)  micron  region). 

Results  of  Task  e) 

With  respect  to  this  task,  no  validation  was  accomplished.  Further, 
the  theory  as  developed  is  not  restricted  to  the  visible  spectrum.  All 
the  assumptions  should  still  be  valid  in  the  near  infrared  region  (0.7 
to  0.9  micron).  The  major  effect  of  wavelength  should  be  a change  in 
certain  parameter  values. 


Measurement  methods  for  determining  atmospheric  Imaging  parameters 


r 


! 


are  discussed  for  the  simple  contrast  reduction  theory  of  Duntley  and 
Middleton  and  for  the  field  theory  development  of  resolution  and  contrast 
dependence  on  all  target  and  atmospheric  parameters. 

Finally,  some  conclusions  about  the  study  are  drawn  and  discussed 
in  terms  of  areas  needing  further  work. 


Rigorous  solution  of  the  general  effects  of  the  atmosphere  on  propa- 
gating electromagnetic  energy  requires  the  solution  of  Maxwell's  equations 
for  a time-varying,  inhomogeneous  medium  subject  to  a given  set  of  boundary 
conditions.  The  atmosphere  is  assumed  to  be  locally  isotropic  since 
measurement  data  tend  to  agree  with  this  assumption.  Even  so,  only 
approximate  solutions  are  available  and  the  validity  of  each  approximation 
method  must  be  checked  against  the  parametric  definition*for  a specific 
application.  For  example,  early  attempts  to  study  nrcpagation  in  a random 


medium  used  a geometric-optics  approximation  which  was  unfortunately 
valid  for  only  very  short  optical  paths.  A later  development  by 
Tatarski  [1]  which  is  based  on  the  Rytov  approximation  and  which  has  been 
widely  used  provided  a much  greater  range  of  validity.  However,  for 
propagation  paths  In  excess  of  one  kilometer,  even  the  Rytov  approxima- 
tion  begins  to  diverge  from  experimentally  derived  propagation  measure- 

» 

ments.  Based  on  the  work  of  Lee  and  Harp  [2],  which  shows  the  Rytov 
approximation  to  be  equivalent  to  the  scatter  of  the  incident  wave  by 
a series  of  random  phase  screens,  the  divergence  from  measured  values  at 
long  optical  paths  is  due  to  a neglect  of  multiple  scattering  effects 
of  small  regions  of  turbulence  along  the  optical  path. 

Of  the  several  methods  developed  for  overcoming  the  limitations 

of  the  Rytov  method,  the  most  promising  is  the  Markov  approximation  [3]- 

[5]  which  has  validity  over  ranges  in  excess  of  hundreds  of  kilometers. 

For  a detailed  comparison  of  the  above  approximations  see  Fante  [6]. 

Three  separate  but  related  phenoma  contribute  to  the  effects  of 

the  atmosphere  on  optical  radiation:  absorption,  scattering  and  refractive 

index  fluctuations.  Absorption  and  scattering  are  typically  related  to 

the  effects  of  atmospheric  gas  and  particle  contituents;  whereas, 

t»  iurbubnce  effect* 

refractive  index  fluctuations  are  typically  relatec^(1.e.  density  gradients, 
temperature,  and  pressure  differences).  Absorption  and  scattering  effects 
on  contrast  between  an  object  and  its  background  is  based  on  work  by 
Koschmieder  [7]  which  was  later  simplified  in  a paper  by  Duntley  [8]  to 
give  the  "two-constant"  theory.  That  Is,  the  reduction  of  contrast  by 
the  atmosphere  is  adequately  represented  by  two  parameters  for  most  "seeing" 


conditions.  The  two  parameters  are  hot  in  fact  constants  but  depend  on 
optical  wavelength  and  on  the  specific  geometry  of  the  observation 
scenario.  Middleton  [9]  present  an  easy  to  read  development  of  the 
Koschmleder  theory  and  of  Duntley's  "two-constant"  theory.  In  addition 
he  presents arguments  for  eliminating  certain  restrictive  assumptions 
made  by  Duntley.  The  work  by  Duntley  and  by  Middleton  still  is  the 
basis  for  contrast  reduction  formulas  cited  in  more  recent  treatises  on 
the  optical  effects  of  the  atmosphere  [10]. 

The  second  body  of  knowledge  found  in  the  literature  treats  the 
effects  of  atmospheric  turbulence  on  MTF  and  typically  makes  no  mention 
of  what  relation  this  has  to  contrast.  Recent  works  [11]-[13],  however, 
include  absorption  and  ambient  light  in  developing  the  EM  theory  approach 
to  imaging  in  an  Inhomogeneous  medium.  The  work  of  Duntley  and  Middleton 
1.n  turn  develops  no  relation  of  contrast  reduction  to  resolution  or  MTF 
effects  of  the  atmosphere.  From  their  definitions  we  know  that  there  Is 
a close  relationship  between  MTF  and  Contrast  Transfer  Function  (CTF). 

So  the  question  which  arises  Is,  "Do  the  two  apparently  different  results 
derived  from  different  approaches  to  the  problem  actually  give  the  same 
Information?"  This  question  will  be  examined  later  in  this  report  after 
a more  detailed  outline  is  given  of  results  of  the  two  different 
approaches. 

The  basic  primer  for  work  relating  to  the  MTF  of  the  atmosphere  Is 
the  book  byTitcirskl  [1],  A later  book  by  the  same  author  [14]  extends 
results  presented  in  the  first  book  and  Includes  a discussion 
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of  the  Markov  approximation  method.  More  recent  works  have  expanded, 
modified  and  given  different  approaches  to  the  basic  work  of  Tatarskl. 

In  most  cases  MTF  of  the  atmosphere  Is  examined  by  first  analyzing  the 
point  spread  function  of  a lens  for  which  the  Input  wave  exhibits  a 
random  log-amplitude  and  phase  variation  across  the  lens  aperture. 

The  MTF  is  then  the  Fourier  transform  of  the  point  spread  function. 

Since  the  log  amplitude  and  phase  are  random  functions,  only  statistically 
averaged  results  are  obtainable.  Additionally  for  short  exposure  times 
the  results  are  such  that  the  effect  due  to  the  atmosphere  Is  not 
separable  from  the  effect  due  to  the  lens.  Fried  [15]  has  presented  a 
development  which  starts  with  fundamental  relations  given  In  Born  and 
Wolf  [16]  and  proceeds  to  final  MTF  using  concepts  and  techniques  which 
are  more  familiar  to  electrical  engineers.  He  develops  expressions  for 
both  short-time  and  long-time  exposure.  Hufnagel  and  Stanley  [17] 
develop  an  optical  transfer  function  for  the  atmosphere  based  on  compu- 
tation of  the  spatial  mutual  coherence  function  for  an  atmospherically 
refracted  wavefront.  Their  result  is  based  on  a long  exposure  time  so 
that  the  MTF  is  accurately  represented  by  an  ensemble  average  with 
respect  to  the  log-amplitude  and  phase  random  functions.  Fried’s  work 
[15]  which  was  done  later  agrees  with  the  Hufnagel  and  Stanley  result 
for  a long-time  exposure  case.  A recent  article  by  Fante  [6]  serves  as 
a good  summary  review  of  the  optical  effects  of  atmospheric  turbulence 
and  offers  a rather  extensive  bibliography. 

An  extension  of  the  monochromatic  theory  for  contrast  reduction  as 
developed  by  Durttley  [8]  and  Middleton  [9]  was  accomplished  and  results 
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are  described  in  this  report.  Example  calculations  show  that  the  error 
incurred  by  using  a monochromatic  representation  for  the  broadband  contrast 
reduction  of  the  atmosphere  Is  very  small  for  typical  imaging  scenarios. 
This  result  is  in  agreement  with  the  opinions  of  other  Investigators  [18] 
that  wavelength  effects  are  minor  effects  over  a relatively  narrow 
spectral  region  (as,  for  example,  the  spectral  region  of  sensitivity  for 
an  imaging  seeker).  A recent  article  [19]  describes  the  effects  of 
the  atmosphere  on  optical  signals  which  exhibit  partial  spatial  coherence 
(such  as  might  be  observed  by  illuminating  a rough  surface  with  a laser); 
however,  the  analysis  is  still  for  monochromatic  radiation. 

Also,  in  this  report,  a detailed  theoretical  development  based  on 
Maxwell's  equations  is  given  which  describes  the  optical  properties  of  ...  . 
the  atmosphere.  The  development  includes  absorption.  From  the  resultant 
atmospheric  Green's  function,  the  propagation  of  complex  spatial  optical 
amplitude  functions  from  the  object  to  the  inage  plane  of  a simple  lens 
is  given.  The  atmosphere- lens  MTF  is  then  computed  from  the  spatial 
Fourier  transform  of  the  optical  intensity  function  in  the  image  plane. 

Contrast  in  the  image  plane  is  computed  in  terms  of  the  definition 
for  universal  contrast  using  the  image  plane  intensity  and  adding  a term 
for  ambient  irradiance.  The  theory  as  developed  gives  a rigorous  treat- 
ment of  the  contrast  transfer  and  MTF  of  an  atmosphere-lens  imaging 
system  for  the  monochromatic  case.  In  terms  of  extending  the  theory  to 
finite  spectral  band  imaging  the  following  apply.  The  development  Is 
general  up  to  a point  at  which  an  assumption  of  monochromaticity  was  made. 
Thus  the  theoretical  development  to  that  point  applies  to  the  broadband 
case.  It  remains  to  extend  the  model  to  give  contrast  amd  MTF  for  the 
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broadband  case. 

Measurement  methods  for  determining  atmospheric  imaging  parameters 
are  discussed  in  Section  5.0  for  the  simple  contrast  reduction  theory  of 
Duntley  and  Middleton  and  for  the  field  theory  development  of  resolution 
and  contrast  dependence  on  all  target  and  atmospheric  parameters. 

Finally,  some  conclusions  about  the  study  are  drawn  and  discussed 
In  terms  of  areas  needing  further  work. 
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II.  CONTRAST  REDUCTION  BASED  ON  A THERMODYNAMIC  APPROACH 

The  development  for  optical  contrast  reduction  by  the  atmosphere  Is 
based  on  work  by  Duntley  [8]  and  Middleton  [s]  and  assumes  the  scattering 
of  radiation  into  and  out  of  a given  optical  path  can  be  adequately 
represented  by  two  constants.  These  two  constants  are  the  atmospheric 
extinction  coefficient,  a,  and  the  ambient  air  light,  The  extinction 
coefficient  includes  scattering  and  absorption  and  is  dependent  on  wave- 
length as  well  as  the  composition  and  size  of  atmospheric  constituents. 
Thus  a will  depend  on  many  separate  mechanisms  which  in  general  are 
difficult  to  predict  theoretically.  With  the  exception  of  the  relatively 
restrictive  cases  of  Rayleigh  and  Mie  scattering  most  data  on  extinction 
coefficient  values  is  based  on  empirical  measurements.  The  meteorologi- 
cal range,  V,  is  defined  as  the  range  at  which  the  contrast  reduction 
due  to  the  atmosphere  is  0.02.  Also  called  the  "visibility  range"  or  just 
"visibility",  V has  units  of  distance  and  has  been  assigned  values  for 
several  subjective  descriptions  of  "seeing  conditions".  For  example  a 
standard  clear  day  is  one  for  which  the  visibility  is  23.5  km. 

A well -presented  and  easily  readable  review  of  the  physics  of  light 
scattering  is  given  by  Harris  [20].  His  definition  for  scattering  is 
"a  process  by  which  radiant  energy  incident  on  an  arbitrary  medium  is 
redistributed  into  one  or  more  angular  directions  of  propagation".  He 
includes  absorption  as  a scattering  mechanism  In  addition  to  refraction, 
reflection  and  diffraction.  Incident  light  will  be  refracted  and/or 


reflected  at  a refractive  index  boundary.  The  angles  and  their  relation- 
ship to  the  refractive  indices  are  given  by  the  Fresnel  fumulae  found  in 
most  elementary  optical  references.  Diffraction  is  the  bending  of 
light  rays  which  pass  close  to  an  edge  or  opening.  For  the  case  of  the 
atmosphere,  light  will  be  diffracted  when  passing  close  to  particles  in 
the  air.  Finally  a number  of  absorption  mechanisms  may  contribute  to  the 
overall  scattering  process.  Resonance  absorption  and  re-emission 
changes  only  the  direction  of  the  radiation;  wavelength  is  unchanged. 
Fluorescence  and  Raman  absorption  on  the  other  hand  result  in  a change 
in  wavelength  for  the  re-radiated  photon.  Total  absorption  of  a photon 
by  a particle  results  in  an  increase  of  the  internal  energy  of  the  particle 
which  essentially  re-distributes  the  radiated  energy  over  a blackbody- 
like spectrum. 

The  rigorous  treatment  of  scattering  is  again  based  on  solutions  to 
Maxwell's  equations.  The  Mie  theory  develops  solutions  for  the  case  of 
small  spherical  particles.  Rayleigh  theory  develops  solutions  for  very 
small  particles  (actually  the  Rayleigh  theory  applies  primarily  to 
molecules  of  the  atmospheric  gases  and  not  to  particles,  and  shows  an 
Inverse  sixth  power  dependence  on  spherical  radius  for  particles  £.1 
micron  for  visible  radiation.)  The  details  of  a scattering  situation  are 
a complex  function  of  particle  size,  particle  shape,  the  complex  index 
of  refraction  and  angle. 

Measurement  of  scattering  cross-sections  from  few  particles  of  a 
given  size  is  extrapolated  to  more  dense  particle  populations  directly 
since  laboratory  and  field  measurements  Indicate  the  validity  of  this 


of  this  approach.  In  other  words,  the  spacing  between  particles  Is  much 
greater  than  the  particle  size  and  the  optical  wavelength  so  that  there 
is  no  Interaction  between  particles.  However,  for  scattering  In  high 
particle  concentrations  such  as  a dense  cloud  this  Is  not  quite  true  and 
the  total  scattering  effect  will  be  smoothed  by  multiple  scattering.  In 
fact  this  effect  Is  more  pronounced  at  optical  wavelengths  than  at  infra- 
red wavelengths  so  that  seeing  conditions  in  a very  hazy  atmosphere  will 
be  better  at  infrared  wavelengths.  The  degree  of  haziness  at  which  this 
advantage  for  the  infrared  becomes  obvious  is 

not  well  defined  and  is  strongly  dependent  on  the  specific  atmospheric 
conditions. 

A.  The  Theory  of  Duntley  and  Middleton 

The  following  development  of  the  contrast  reduction  effects  of  the 
atmosphere  is  based  on  the  original  work  of  Duntley  as  later  modified  by 
Middleton.  The  purpose  here  is  to  g1*'e  a brief  theoretical  development, 
identify  the  general  result  and  discuss  special  cases,  in  particular  the 
sky-to-ground  observer. 

The  model  for  contrast  reduction  due  to  the  atmosphere  for  a sky-to- 
ground  observer  is  shown  in  Figure  1.  For  a differential  lamina  of  air 
perpendicular  to  the  path  of  sight  we  have 

L(r)  + dr  = L(r  + dr)  (1) 

where:  L(r  + dr)  is  the  radiance  out  of  the  lamina  toward  the  observer 
and  dL(r)/dr  is  the  change  in  path  radiance  in  the  lamina. 

Clearly,  radiance  can  be  lost  or  gained  in  the  lamina  dr,  even  In 
the  absence  of  sources  within  dr  which  is  the  case  assumed  for  the  atmosphere. 


. ..... *•*  ^ ~ 
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Then  dL(r)/dr  depends  on  two  things:  1)  radiance  lost  through  extinction 
(scattering  and  absorption)  and  2)  radiance  gained  as  air  light  from  the 
surrounding  medium  which  has  been  scattered  into  the  path  of  sight. 


Figure  1.  Model  for  contrast  reduction. 


Then 

= -0(r)L(r)  + La(r)  (2) 

where:  a = extinction  coefficient 

a “ a + 8 

a ■ absorption  coefficient 
6 ■ scattering  coefficient 

la(r)  Is  the  air  light  scattered  Into  the  path  of  sight 
In  arriving  at  the  differential  equation  In  (2)  a number  of  assumptions 
are  Inherent.  First  It  Is  assumed  that  the  light  scattered  by  a lamina  Is 
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directly  proportional  to  the  incident  light.  This  assumption  Is  implied 
by  the  argument  that  there  is  no  Interaction  among  molecules  or  particles 
In  terms  of  scattering.  Certainly  for  a reasonable  range  of  particle 
densltltes  and  incident  radiation  levels  this  superposition  theory  Is 
expected  to  hold  true.  Other  assumptions  in  the  development  of  (2)  are 
that  the  object  be  small  in  comparison  to  the  range  R and  that  second  order 
scattering  of  light  from  the  object  into  the  path  of  sight  is  inseparable 
from  and  a part  of  the  term  for  air  light,  L_ ( r) . Additionally  It  should 
be  emphasized  that  the  terms  in  (2)  are  dependent  not  only  on  r as  Indicated 
but  on  optical  wavelength  as  well.  Without  knowing  the  specific  functional 
form  of  the  "r" -dependence  we  can  go  no  further.  Both  Duntley  and  Middleton 
provide  a convenient  way  for  avoiding  this  difficulty.  It  is  based  on  the 
assumption  that  o(r)  and  La(r)  both  have  the  same  functional  dependence 
on  r so  that 


o(r)  - o0  f (r)  (J) 

L*(r)  = La(0)f(r) 

where  f(r)  is  some  unspecified  function  of  r. 

Oq  and  La(0)  are  constants  and  represent  their  respective  values  at 
r = 0.  With  this  assumption  (2)  becomes 


where  the  right  hand  side  of  (4)  is  designated  R and  is  called  the  "optical 


RADIANCE  EFFECTS  OF  THE  ATMOSPHERE 


t-Q  Is  the  Inherent  object  radiance 

Lp  is  the  apparent  object  radiance  at  slant  range  R.  Note 
that  Lr  has  two  terms.  The  first  term  represents  an  increase  in  apparent 
radiance  due  to  air  light  being  scattered  into  the  optical  path  and  the 
second  term  shows  how  the  inherent  object  radiance  is  reduced  by 
atmospheric  extinction.  Equation  (5)  is  the  basic  result  of  Middleton  and 
Duntley.  Other  expressions  which  relate  contrast  of  the  object  with  a 
background  and  for  various  viewing  geometries  occur  more  often  in  the 
literature  than  (5). 

For  example  if  we  consider. a localized  background  adjacent  to  the 
object  (remember  that  the  initial  assumptions  required  the  object,  or 
background,  to  be  small  compared  to  R)  then  its  apparent  radiance  at  R 
is  given  by: 


And  if  we  define  contrast*  as 


C $ 

L 

then  the  inherent  contrast  is 


and  the  apparent  contrast  at  R is 


(7) 

(8) 

(9) 


♦There  are  several  variations  on  the  definition  for  contrast.  See 
Appendix  B. 
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Substituting  (5)  and  (6)  Into  (9)  gives 


CONTRAST  TRANSFER  OF  THE  ATMOSPHERE 


Equation  (11)  is  the  general  result  for  contrast  reduction  by  the  atmo- 
sphere given  by  Middleton  and  by  Duntley.  This  result  was  developed  for 
the  air-to-ground  observer  case;  hpwever.  It  Is  valid  for  the  ground-to- 
air  case  except  that  the  air-light  contribution  in  equation  (6)  should 
be  different  for  the  two  viewing  directions. 

For  certain  special  cases  where  the  background  radiance  is  independent 
of  range  (the  commonly  cited  example  is  for  horizontal  viewing  against 
a sky  background)  we  have  that  Lq'  = LR'  and  (11)  is  written  as 


SPECIAL  CASE:  RANGE  INDEPENDENT  BACKGROUND  RADIANCE 
An  alternate  expression  for  the  air-to-ground  case  is  given  by 


substituting  (6)  into  (11)  for  LR*  to  give 


c« 

*0  ■ 


- (13) 


mm 

which  is  further  modified  by  recognizing  ^21  , £hc 


reflectivity  of  the  background  (see  Middleton  [9]  p-71  for  a discussion). 
Thus 


Is  a valid  expression  which  can  also  be  shown  to  be  equal  to 


(14) 
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L 


H’  - ¥ H <,5> 

-o0R 

where  T * e is  the  atmospheric  transmittance  over  optical  path  IT. 
B.  Extension  to  a Finite  Spectral  Band 


In  evaluating  the  usefulness  of  theoretical  expressions  for  contrast 
reduction  by  the  atmosphere  in  terms  of  applicability  for  determining 
performance  of  imaging  seekers  or  other  broadband  optical  systems,  it 
becomes  necessary  to  examine  the  effert  of  spectral  dependence.  Some 
method  for  reconciling  a monochromatic  theory  with  a wideband,  spectrally 
non-trivial  application  such  as  Is  true  for  visible  and  near-infrared 
imaging  seekers  must  be  found. 

From  the  work  of  Duntley  and  Middleton,  the  spectral  radiance 
transfer  characteristic  of  the  atmosphere  is  given  for  monochromatic 
radiation  by 


L „ <R) 


^JoT-  [1  - TJ  * La(0}Tx 


(16) 


where  L^(R)  Is  the  spectral  radiance  at  slant  range  R due  to  a source 
radiance  LA(0)  at  zero,  oA(0)  is  the  atmospheric  extinction  coefficient 
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at  zero  (including  scattering  and  absorption)  for  a given  wavelength, 
La>x(0)  is  the  spectral  radiance  of  the  air  light  and  Is  the  atmospheric 
spectral  transmittance  given  by 


(22) 


where  L^x(0)  = J lx(0)dx.  From  (20)  and  (21)  we  get 
AX 


cax<r> 


K«»l 

L^J 


'AX 


Equation  (22)  is  in  the  same  form  as  the  general  result  derived  by 
Duntley  for  the  monochromatic  case,  where  TAX  is  now  the  mean  broadband 
atmospheric  transmittance.  We  can  equate  T^x  with  exp[-oAX(0)EAX] 
where  oAX (0)  can  be  called  the  mean  broadband  extinction  coefficient  of 
the  optical  path.  However,  o^x (o)  is  not  a simple  average  value  within 
the  spectral  band  of  interest  but  is  given  by 

f {[Lx(0)  - LJ(0)]T,}dX 

Ax 


°AX(°) 


A 

R 


In 


AX 


L 


(Lx(0)  - L'(0))dX 


(23) 


L j AX 

The  result  in  (23)  is  easily  evaluated  by  computer  if  all  the  spectral 
data  are  available.  For  the ‘special  case  of  flat  spectral  emission  by 
the  object  and  background  we  get  the  simpler  result  which  is  Independent 
of  source  and  background  radiance  levels 


AX 


(0) 


AX 


exp[-o, (0)Rx]dx 


(24) 


Further,  if  we  look  at  the  derivation  for  R in  Duntley* s paper  [8] 
we  have  the  following 


and 


o.(r)  * o,(0)f(r) 


r R 

R * / f(r)dr 
J n 


(25) 
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where  the  unspecified  function  f(r)  was  assumed  to  give  range  dependence 
for  both  the  extinction  coefficient  and  the  air  light,  Lfl  x(r).  It 
seems  reasonable  to  argue  that  the  spectral  dependence  for  ox(r)  is 
contained  in  the  ax(0)  term  so  that  f(r)  and,  hence  R,  are  not  functions 
of  wavelength.  With  this  argument  the  x subscript  can  be  dropped  from 
R so  that 


\ ■ R 


Then  (23)  becomes 


°ax(°)  s ln 


/ (L  (0)  - LM0))dx 
/ax  _ 

/ t[Lx(0)  - L'(0)]Tx}dX 


or  for  the  case  of  flat  target  and  background  spectral  radiance  (24)  is 


given  by 


°AX(0) 


■ » ,n  ( L 


exp[-ax(0)R]dx 


The  result  in  (28)  is  compared  with  the  value  of  extinction  coefficient 
obtained  as  a simple  average  value  across  AX  by  the  ratio 


c . 

oautTO t 


"s'"  ^ -Cx  «Pt-«x(0)R]dx| 


The  comparison  is  made  by  choosing  a model  for  ox(0)  in  the  visible 
spectrum.  Several  authors  (for  example  [21]  and  [22]  have  suggested  models 


of  the  form 


o,(0) 


(i)" 


where  oq.55^0^  *5  the  extinction  coefficient  at  a wavelength  of  0.55 
micrometers  for  a given  visibility,  "a"  Is  0.55  or  0.61  micrometers  and 
"q"  is  a constant  in  the  range  0.7  to  about  3.0  for  typical  visibility 
values. 

The  effect  on  atmospheric  transmittance  of  a broadband  theory  Is 
further  compared  with  monochromatic  theory  using  average  values  for 
ax(0)  by  the  ratio 


TAX  = exp[-o&x(°)R] 
*^6  exp[-oAV6(0)R] 


(31) 


Values  for  p are  shown  in  Figure  2 for  various  visibilities  as  a 
function  of  slant  range.  Parameter  values  a = 0.55  and  q * 0.7  were 
used  in  the  computation. 


Figure  2.  Ratio  of  broadband  to  average  transmittance  for 
various  optical  slant  ranges  and  visibilities. 
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The  result  in  (22)  is  dependent  on  an  equivalent  broadband  extinction 


i 


coefficient  o&x(0),  which  is  dependent  on  the  spectral  radiance  levels 
of  the  target  and  background  as  well  as  the  spectral  extinction 
coefficient  <>x(0).  Rigorous  computation  of  the  broadband  contrast 
reduction  is  easily  achieved  by  computer  if  all  the  spectral  data  are 
available. 

A special  case  where  target  and  background  radiances  are  spectrally 
flat  yields  the  results  in  Figure  2.  For  high  visibility,  little  error 
Is  encountered  in  the  atmospheric  transmittance  by  using  an  average  value 
for  the  extinction  coefficient.  For  poor  visibility  at  longer  ranges 
the  effect  becomes  more  pronounced.  It  is  emphasized  that  these  results 
are  for  simplified  (spectrally  smooth)  models  for  target/background 
radiances  and  spectral  extinction  coefficients.  Models  which  included 
distinct  spectral  characteristics  should  have  a greater  effect  on  the 
approximation  error,  p. 

A major  assumption  made  in  the  development  Is  one  made  by  Duntley 
and  by  Middleton  in  assuming  that  both  the  extinction  coefficient  ax (R) 
and  the  ambient  light  scattered  into  the  optical  path  L . (R)  have  the 
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same  functional  dependence  on  range  R.  Extinction  coefficient  includes 
effects  of  both  scattering  and  absorption  and  depends  on  wavelength, 
particle  size  (r * ) and  refractive  index  (m).  Additionally  the  scattered 
and  absorbed  radiation  exhibits  an  angular  distribution  which  is  dependent 
on  these  same  parameters.  Examples  of  the  dependence  of  scattering  and 
absorption  on  the  size  parameter  a » 2*r'/X  for  various  values  of  m and 
scattering  angles  are  given  in  the  excellent  sumnary  paper  by  Harris  [20]. 


The  "air  light"  term,  L . (R),  represents  radiation  scattered  into 
the  line  of  sight  by  particles  along  the  path.  This  term  is  then 
dependent  primarily  on  the  scattering  properties  along  the  path  plus  the 
angular  distribution  of  ambient  radiation.  It  is  easy  to  visualize 
conditions  for  which  o(R)  and  L,  , (R)  have  significantly  different  range 
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dependencies. 

However,  the  major  thrust  of  this  development  was  to  demonstrate 
that  for  certain  imaging  applications  where  sensitivity  over  the  visible 
spectrum  exists,  some  mechanism  for  including  spectral  dependence  in 
contrast  reduction  calculations  is  desirable. 

In  summary,  the  following  points  are  made  with  respect  to  Figure  2. 

a)  A value  of  TAX/T^yg  greater  than  unity  shows  that  the  atmospheric 
transmittance  is  larger  using  the  broadband  model  than  if  an  average 
extinction  coefficient  is  used.  This  improves  contrast  slightly. 

b)  The  method  for  finding  T^yg  is  more  rigorous  than  some  averaging 
techniques,  eg.,  one  method  Is  to  use  the  value  of  a at  the  central  wave- 
length in  the  spectral  passband. 

c)  The  dependence  of  o on  wavelenth  was  modeled  as  a smooth  function'. 
In  some  spectral  regions,  o exhibits  more  distinct  wavelength  dependence. 

d)  For  the  case  examined.  Figure  2 shows  that  T^yg  and  TAJl  are 
essentially  equal  for  slant  ranges  less  than  the  visibility  range. 

e)  The  results  given  are  for  an  assumed  flat  spectral  radiance  for 
both  the  target  and  background. 

f)  Any  changes  In  b.),  c.)  and  e.)  above  will  affect  the  result  shown 
in  Figure  2 . The  exact  effect  is  strongly  case  dependent. 
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III.  OPTICS  OF  THE  ATMOSPHERE  BASED  OH  ELECTROMAGNETIC  FIELD  THEORY 
In  this  section,  beginning  with  Maxwell's  equations,  a detailed 
development  of  the  theory  behind  effects  of  the  atmosphere  on  optical 
signals  is  given. 

First,  in  3.1,  the  vector  Helmholtz  wave  equation  describing  vector 
electric  and  magnetic  fields  for  the  optical  signal  is  derived  from 
Maxwell's  equations,  several  vector  Identities,  assumptions  as  to  the 
nature  of  the  optical  path  and  from  Fourier  theory  for  a monochromatic 
optical  waveform.  Next  the  scalar  monochromatic  theory  is  derived  from 
the  vector  theory  which  gives  the  scalar  Helmholtz  wave  equation.  It 
describes  the  spatial  behavior  of  an  optical  complex  amplitude  propagating 
in  an  isotropic,  linear,  slightly  inhomogeneous,  slightly  absorbing 
medium  (l.e.  the  atmosphere).  Solutions  to  the  scalar  Helmholtz  equation 
subject  to  a set  of  boundary  conditions  are  then  approximated.  A number 
of  methods  have  been  used  to  approximate  a solution  to  the  Helmholtz 
equation.  Several  of  these  are  outlined  briefly.  A conceptually  satis- 
fying solution  is  found  using  a Green's  function  approach.  The  solution 
is  the  Rayleigh-Sommerfeld  integral  extended  to  included  inhomogene ites 
and  absorption.  A further  assumption  produces  the  modified  form  of  the 
Rayleigh-Soninerfeld  Integral  or  the  extended  Huygen-Fresnel  Integral. 

This  method  Is  given  in  considerable  detail . 

A.  Maxwell's  Equation*  and  the  Inhomogeneous  Scalar  Helmholtz  Equation 
Using  rationalized  MKS  units,  we  start  with  Maxwell's  equations  [23] 
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7 x H * D + T 

7 x E * -t 

V • D ■ p 

7 • B ■ 0 


Isotropic,  Linear  Medium 


U = c(x,  y,  z)F 
B = y (x , y,  z)H 
J = o(x,  y,  z)F 


divide  (33)  by  y and  use  (37) 


- 7 x E = -H 
v 


Take  curl  of  (39) 


v x (^  v x E)  = -7  x H 


differentiate  (32)  wrt  time 


7 x H = D + J 


use  (36)  and  (38)  in  (41)  and  if  e and  o not  fens,  of  time* 


7 X H * cE  + oE 


where  x,  y,  z dependence  of  e,  a and  y will  be  left  implicit  for 
convenience. 


+ In  general,  e and  o,  as  well  as  the  refractive  index  to  be  defined 
from  them,  will  depend  on  x,  y,  z and  on  "t".  However,  this  time 
dependence  is  usually  assumed  to  be  associated  with  atmospheric  winds 
and  according  to  "Taylor's  frozen  flow  hypothesis"  is  statistically 
dependent  on  [r  - V(r)t]  where  r ■ (x,  y,  z)  and  V(r)  is  the  local  wind 
velocity.  For  the  development  here,  we  take  the  approach  of  most  authors 
and  ignore  the  time  dependence  in  obtaining  approximate  solutions  to  the 
wave  equation.  Time  dependence  is  considered  later  in  the  development  in 
terms  of  specific  exposure  times  and  expected  effects  on  MTF. 
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Substitute  (42)  into  (40) 

V X (i  V X f ) * -(eF  + oE)  (43) 

identity 

7 x (uV)  » uv  x V + vu  x V (44) 

u is  a scalar 

V is  a vector 

Using  (44)  in  (43)  where  u = 1 and  V =•  7 x E 

-7  x 7 x E + 7(1)  x V x E = -(cE  + oE)  (45) 

V w 

identity  7 x 7 x V = 7 (7  • V)  ■ v2V  (46) 

using  (46)  in  (45)  gives 

1 7(7  • E)  - 1 V2E  + 7(1)  x 7 x E = -(eE  + oE)  (47) 
From  (34),  (36)  and  using  the 

identity  div  uV  = 7 • (uV)  = u7  • V + V • 7u  (48) 

7 • D * 7 • (eE)  = e7  • E + E • 7e  = p (49) 

Solving  (49)  for  7 • E 

7 • E = 1 [p  - E • 7e]  (50) 

Substitute  (so)  into  (47) 

77  (l(p  - E • 7e) ) - ^72E  + 7(1)  x 7 x E = -(eE  + oE)  (51) 
Multiply  (51)  by  p 

7{1  (p  - E • 7e]}  - 72E  + p 7(1)  X 7 X E + u [eE  + oE]  = 0 (52) 

or  rearranging  and  changing  sign 

72E  - peE  - poE  - p7(l)  X 7 X E + 7 [If  . 7e]  - 7[^\)  = 0 (53) 
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Result  in  (53)  is  correct  for  a linear  isotropic  medium  (i.e.  all  terms 
in  E are  of  power  one,  p & p & e & a may  depend  on  x,  y,  & z but  will 
be  same  at  a given  point  for  all  components  of  the  vector  F. 

An  alternate  form  of  (53)  is  achieved  by  recognizing  that 


y(p>  - * J2  ^ 

y 


v(logp)  = Vp  = ^7p 

op  P 

From  (54)  and  (55)  we  see  the  term  in  (S3) 


[-P  v (1)3  becomes  -p(-  -^Vp)  = -j7p  = v(log/$  (56) 

p 

Further  since 

* * ^ = J A • vu  (57) 

and  v(log  e)  = (58) 

C 

the  term  1 E • 7c  becomes  E • 7 (log  e)  (59) 

Using  (56)  and  (59)  in  (53)  gives 

Vector  D.E.  for  electric  field  T,  Linear,  Isotropic  Medium 

2—  *•_  • __ 

7 E - peE  - paE  + V ( 1 og  p)  X 7 X E + 7[E  • 7(log  e)]  - v(£)  = 0 


SPECIATCASE  - VACUUM  (60) 

p * constant  , e * constant  -►  7(log  p)  = 7(log  e)  * 0 


o = 0 


. p * 0 


zero  conductivity 


source- free 


Result  is  Vector  Helmholtz  equation 

2 •• 

7 r - pcF  = o 


I 


SPECIAL  CASE  - ATMOSPHERE 

Conditions  (1)  Weakly  Inhomogeneous  - l.e.  w,  e,  A o show  slight 
variations  with  x,  y & z;  however,  the  variation  path  length  Is  large 
compared  to  spatial  variation  in  F (>v  wavelength) 

! Result  of  (1):  In  (60)  v(log  v)  A v (log  c)  are  small  and  approxi- 
i mated  as  zero. 

(2)  Weakly  absorbing  (conducting)  - i.e.  charge  density  p induced 
'by  absorption  of  energy  is  small  and  varies  slowly  with  x,  y,  z. 

Result  of  (2):  v(^)  is  small  and  approximated  as  zero. 

Then  eqn  (60)  becomes 


n . ••  « 

v E - peF  - poF  * 0 

Weakly  in-homogeneous,  weakly  absorbing 
medium  (linear  and  isotropic  also) 

(62) 


To  handle  time  derivatives,  take  temporal  Fourier  transform  of  E(x,y,z,t) 
E(x,y,z,v)  , v is  optical  frequency  (Hz) 
l.e.  the  defining  equations  for  temporal  transform  are 


Temporal 

Fourier 

Transform 

Pair 


- . _ j2irvt 

E(x,y, E(x,y,z,t)e  dt 

— A /•»  * -j2irvt 

E(x,y,z,t)— E(x,y,z,v)e  dv 


(63) 

(64) 


If  (64)  is  differentiated  wrt  time  once  & twice,  the  result  is 

• • A - j2wvt 

£ s / C-j2nv  E(x,y,z,v)]e  dv  (65) 

— 00 

•i  o-  -j2nvt 

E * / C(-j2»vrE(x,y,z,v)]e  dv  (66) 


I 


I 


Thus  we  have  the  following  transform  pairs 
E(x,y,z,v)  <»  E(x,y,z,t) 

A ^ 

-j2irv  E(x,y,z,v)  E(x,y,z,t) 

(-j2nv)2  Elx,y,z,v)  » F(x,y,z,t) 

Using  (67)  and  taking  the  Fourier  transform  of  (62)  gives 

V2E  + yc(2irv)^  E + pa(j2uv)E  = 0 


[ * ♦ (2irv)S(e  + j^-)]E 


E=E(x 


-Vector  Helmholtz  Wave  Equation- 
Weakly  inhomogeneous 
Weakly  Absorbing 
Isotropic 
Linear  Medium 


The  term  in  parenthesis 


*e  + j 2*J  = G 


can  be  thought  of  as  a complex  permittivity  in  which  case  (68)  Is  identical 
to  the  temporal  Fourier  transform  of  the  vacuum  equation  (61). 

NOTE:  In  general  for  a non-dispersive  medium  e does  not  depend  on  optical 

A 

frequency  v;  however,  as  seen  by  (69)  the  complex  permittivity  e does 
depend  on  v.  This  is  another  way  of  saying  that  optical  absorption 
(represented  by  the  imaginary  part  of  e)  i$  wavelength  dependent  which 
has  been  verified  by  experiment. 

With  (69),  (68)  becomes 


[v2  + (2wv)2  pe]  E = 0 

Vector  Helmholtz  eqn  for  practical  atmosphere 


Solutions  are  sought  to  (70) 


First  it  is  recognized  that  (70)  describes  only  the  electric  field 
behavior  of  an  electromagnetic  wave.  It  can  be  shown  by  a development 
similar  to  that  above  that  the  magnetic  field  vector  H is  described  by 
a similar  second  order  D.E. 

Second,  note  is  made  that  both  H and  E can  have  three  scalar  compo 
nents,  eg.  in  rectangular  coordinates 


Ex(x,y,z,v)  Hx(x,y,z,v) 
Ey(x,y,z,v)  Hy(x,y,z,v) 
E2(x,y,z,v)  Hz(x,y,z,v) 


(71) 


must  all  be  solved  simultaneously  subject  to  a set  of  boundary  conditions. 
Each  component  in  (71)  is  described  by  an  equation  of  the  form  of  (70). 
Finally,  simultaneous  solution  of  the  components  in  (71)  has  been  achieved 
for  a very  few  special  cases  of  little  practical  interest.  For  cases  of 
practical  interest  we  develop  the  following  scalar  theory. 

Postulate  a scalar  field  $(x,y,z,t)  and  its  temporal  Fourier  trans- 
form $(x,y,z,v)  such  that  $ "represents"  a real  optical  field.  $ may  be 
thought  of  as  the  scalar  amplitude  of  a vector  field  such  as 


(72) 


Thus  letting  <t>  represent  the  optical  field  we  find  that  $ satisfies 
equation  (70) 


Practical 

Atmo- 

sphere 


[v2  + (2irv)2  lie]  $(x,y,Z,v)  = 0 


Scalar  Helmholtz  wave  equation, weakly 
inhomogeneous,  weakly  absorbing, linear, 
isotropic  medium 


(73) 
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Definition  of  some  terms 


k ^ j = wavenumber  (m* 1 


) 


v * velocity  (m/sec) 


c = speed  of  light  (m/sec) 

1_ 
je 

v = 1 / -/  s“ 

\ue 

v = frequency  (sec*1)  = y = kc 

n = refractive  index  ^ - = CV^" 

„ c 

n = complex  refractive  index  = * = 

Then  (2irv)2w2  = (2ir)2k2n2 
Using  (74)  in  (73)  gives 


[v2  + (2*)2k2n2]  «(x,y,z,v)  = 0 


(74) 


(75) 


B. 


Alternate  form  of  (73) 

Solution  Methods  for  the  Helmholtz  Equation 

Solutions  are  sought  to  equation  (39).  Consider  first  the  mono- 
chromatic case,  i.e.,  a real  monochromatic  scalar  field  is  represented 
by 

. -j2nv  t . J2itv  t 

*(x,y,z,t)  = \j>(x,y,z)e  0 + **(x,y,z)e  0 (76) 

whose  transform  is  (with  respect  to  time) 


*(x,y,z,v0)  = ^(x,y,z)6(v-v  ) + **(x.y,z)6(v+v  ) 


(77) 


It  is  verified  that  41  given  by  (77)  satisfies  the  scalar  Helmholtz 

wave  eqn.  (75)  if  ij<(x,y,z)  satisfies  the  scalar  Helmholtz  wave  equation,  i.e.. 


[V2  ♦ (2n)2k2n2]  i(x,y,z)  = 0 


(78) 


r* 
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where:  *(x,y,z)  is  the  complex  spatial  amplitude  function  for  the, mono- 
chromatic scalar  wave  *(x,y,z,v# ). 

Thus  for  the  monochromatic  case,  solution  to  (75)  is  accomplished 
by  finding  a solution  to  (78)  for  4>(x,y,z),  then  substituting  into  (76). 
Approximation  methods  are  found  in  the  literature  which  typically 
characterize  the  index  of  refract i on  (n)  as  follows 


n(x,y,z)  = nQ  + n^x.y.z)  (79) 

where  nQ  » n^ 

nQ  is  the  turbulence  free  refractive  index  of  the  atmosphere  and  n.j  is 
the  fluctuating  component  of  n.  Further  it  is  assumed  that  n^(x,y,z)  is 
a zero-mean,  gaussian  random  variable  so  that  it  is  completely  characterized 
by  its  second  moment,  or  autocorrelation  function. 

1.  Method  of  small  perturbations 

From  equations  (78)  and  (79)  where  ng  is  set  equal  to  unity  and  using 
kQ=  2 irk  we  have 

v2*  + ko2(l  + n-j  )2i|/  = 0 (30) 

A solution  <i>  = ij/g  + ij>i  is  sought  where  ifig  is  the  unperturbed  portion 

2 2 

of  the  wave  solution  and  satisfies  \r<j/0  + kQ  = 0.  Substituting  * = *0  + *i 
into  (80)  and  clearing  gives 

+ ko*l  + 2niko(*0  + *1>  + konl2(*o  + *1>  = 0 (81) 

2 

The  last  term  in  (81)  is  of  order  n-j  and  is  neglected;  further  the  method 
of  small  perturbations  assumes 


1*1 1 « l*0l  (82) 


T 


so  that  equation  (81)  is  written  as 

A + k02<h  * -2"ik02*o 

Solutions  to  the  unperturbed  wave  equation  have  the  form 


*o“V 

where  AQ  and  SQ  are  the  amplitude  and  phase. 

Then  if  we  express  the  solution  + <pi  in  the  form 

* “ *o  + *1  * AelS 

then  Tatarski  shows  that  for  |^|  « l^j  we  have 


log  ip  = log  A + iS  = log  AQ  + iSQ  + ^ 
quating  real  and  imaginary  terms  we  get 


S - Sq  = Im 


Thus  we  have 


*1  = e 


X + iS1) 


= log  x e 


where: 


III * 


1k0|r-r’| 
)»0(r‘)e  0 


Is  the  small  perturbation  solution  to  (83)  In  terms  of  n^  and  <|>0,  where 
r • (x,y,z). 


♦g(r)  which  must  be  known  to  solve  equation  (86)  Is  easily  found  for 
plane  monochromatic  waves,  l.e.  for  a wave  propagating  In  the  z direction 


*0(r)  * V 


1kQz 


(87) 


and  (86)  becomes 


k0K  fff  -Hc0(z-z*) 

>itr)  ■ -kj]J^r">* 


" lr-r'1 


dr' 


(88) 


Further  simplifications  to  (88)  are  made  based  on  small  scattering  angle 
assumptions. 

2.  The  Rytov  method 


This  method  uses  a different  approach  to  solving  (80)  which.  Instead 


of  requiring  smallness  of  compared  to  Vq,  uses  the  less  restrictive 
condition  that  spatial  fluctuations  In  phase  and  amplitude  due  to 
refractive  Index  fluctuations  be  slowly  varying  with  respect  to  the 
optical  wavelength;  that  Is, 


|vSj  | « 2ir 

|vx|  « 1 

The  approach  is  to  find  a solution  * to  (80)  where 
♦ - e X ♦ IS  . eB(r) 
where  B(r)  *■  BQ(r)  + Bj(r) 

and  IvBjl  ««  |vb0|. 


(89) 


(90) 
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It  can  be  shown  by  substituting  (90)  into  (80)  and  e into  (80)  that 
Bq  satisfies  the  vacuum  equation 

(VB0)2  ♦ k02  - 0 (91) 

and  that  B^  satisfies  approximately 

v2Bj  + 2vBg  • vB|  + 2k02n1  * 0 (92) 

The  solution  to  (92)  is  given  by  the  form  in  (86)  so  that  the  method  of 
small  perturbations  and  the  Rytov  method  give  the  same  results.  However, 
The  Rytov  method  has  a wider  range  of  validity.  Hufnagel  and  Stanley 
point  out  that  even  the  Rytov  method  is  based  on  an  assumption  which 
fails  for  |B|  » 1.  Further,  they  develop  a solution  based  on  a statistical 
approach  which  gives  the  mutual  coherence  function  without  finding  an 

explicit  solution  for 

3.  The  Harfcov-PT«bollc  approximation 

Starting  with  an  optical  signal  propagating  In  the  Z direction 
i(kftz-wt) 

E * ♦(x,y,z)e  0 (93) 

the  wave  equation  (80)  is  then  approximated  by 


^ * $ * *,ko  H ♦ *koV  ■ 0 <"> 

This  is  the  so-called  parabolic  approximation  to  the  wave  equation  and 
is  applicable  if  the  beam  spread  angle  for  the  propagating  wave  is  small. 
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In  other  words  with  respect  to  the  succeeding  figure  we  have 

|r-r'|  = y(x-x')Z  + (y-y')Z  ♦ 
where  | x-x* | « z* 

ly-y'l  « Z ' 


Figure  3.  Small-angle  scattering. 


The  solution  to  (94)  can  be  written  as  [13] 


<l«(x,y,z)  = 


fa  f /dx'dy'  *0(x\y',0) 

(96) 

[(x-x1)2  + (y-y')2]  + u(x,x'  ,y,y')} 


where  ^(x*  ,y' ,0)  is  the  field  in  the  z's0  plane,  4>(x ,y ,z ) is  the  field 
in  the  z'=z  plane,  and  u(x,x',y,y')  is  the  random  part  of  the  complex 
phase  of  a spherical  wave  propagating  from  point  (x'.y'.O)  to  point 
(x,y,z).  Equation  (96)  without  the  dependence  on  u(x,x',  y,y')  is  simply 
the  Huygen- Fresnel  formula  which  is  an  approximate  solution  to  the  scalar 
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Helmholtz  wave  equation  in  a homogeneous  medium.  Thus  the  result  in 
(96)  can  be  considered  to  be  an  extension  of  the  Huygen-Fresnel  formula 
to  include  media  with  small-scale  random  fluctuations.  Additionally, 
(96)  is  conveniently  expressed  in  terms  of  linear  systems  theory  since 
it  relates  fields  in  the  input  plane  (x'.y'.z^O)  and  in  the  output 
plane  (x.y.z),  i.e. 


t.y.z)  = JJ  dx'dy' 


*0(x\y',0)  h(x,x',y,y') 


1 K 

where  h(x,x*,y,y*)  = -V  e.{  7T  C(x’x' + (y^' )23  + u(x,x*  ,y,y* )]} 

Z*1z  (98) 

Shapiro  [24]  has  shown  that  h(x,x',y,y')  is  reciprocal,  that  is,  the 

input  and  output  planes  can  be  reversed  and  (97)  still  holds.  Using 

Harkov  approximation  methods,  the  moments  of  i/i(x,y,z)  can  be  found  to 

yield  the  spatial  coherence  function  of  the  propagating  field. 

M(p1,p2,z)  = 0(p-|,z)  <|<*(p2»z)  > (99) 

where  M is  the  spatial  coherence  function  at  z'=z,  the  < > brackets 
denote  an  ensemble  average  and  p^  = (x^.y.).  Using  (97)  in  (99)  gives 


N(pj  »P2  = /A  i*  d*i  JJ dx2'dy2'  <|»q(Xj 1 ,y]  1 ,0) 


(100) 


< h(x1,y1,x1',y1,)h*(x2,y2,x2,,y2,)> 

where  for  the  homogeneous  medium,  < h(p^  tPj' )h*(p2,p2' ) > is  easily 
shown  to  be 

2 

exp{-^  [(p1-p1’)2  - (p2-p2')2]  } (101) 


: - y~ 


and  for  the  fluctuating  medium  Is  shown  not  so  easily  [6]  to  be 

< h(pj,pj')h*(p2.p2')  > exp  { "2J  C (p-j-p-j  ' )2  " (P2"P2*)2^ 

• T2-  / *' H [*••  4 •-.*•>  0 - 9] 


where  H(z,€)  « 1.88  C_z(i)|«|- 


[,-°.«,5(ifi)V3 


002) 

(103) 


Cn  Is  the  refractive  Index  structure  coefficient 
Lq  Is  a measure  of  the  largest  distance  over  which  fluctuations 
In  the  refractive  index  are  still  considered  to  be  correlated. 
Another  way  of  looking  at  what  equation  (100)  means  is  given  by 
recognizing  that  the  function 


M(p1'.p2,,0)  = <h.(p1,,0)  **(p2\0)  > 


(104) 


Is  a measure  of  the  spatial  coherence  in  the  Input  plane  and  equation 
(100)  then  gives  the  propagation  of  this  coherence  function  through  a 
medium  described  by  h(x,y,x',y')  which  is  the  point  response  for  the 
medium.  That  h(x,y,x',y')  gives  the  response  at  (x,y,z)  to  a point 
source  at  (x'.y'.O).  The  relationship  of  spatial  coherence  to  resolution 
or  MTF  is  typically  derived  in  terms  of  the  Fourier  tranform  of  the 
intensity  distribution  in  the  focal  plane  of  a simple  lens  which  is 
illuminated  by  an  atmospherically  distorted  plane  wavefront.  The 
problem  with  this  approach  is  that  for  the  short-exposure  case,  effects 
of  the  atmosphere  are  not  separable  from  effects  of  the  lens.  A 
development  by  Fried  [15J  is  given  in  Section  4 for  MTF  effects  of 
the  atmosphere  for  both  long  and  short-term  exposure. 


4.  The  Hurgcn-Frcgnel  method 

The  Green's  function  method  [25]  is  used  to  find  a solution  to  (78). 
Consider  a volume  V enclosed  by  surface  S.  If  * has  continuous  1st  and 
2nd  order  spatial  partial  derivatives  within  V and  on  S;  further,  if  G 
is  any  other  function  which  satisfies  the  same  continuity  conditions  on 


then  Green's  theorem  yields 


fffy  [iv2G  - GV2 


*]  dV  *ffs  [*  §£  - G||]ds 


(105) 


GREEN'S  THEOREM 


n is  derivative  normal  to  S in  outward  direction.  Now  for  a propagating 
optical  field,  we  are  interested  in  solving  the  boundary- value  problem, 
given  the  field  4>(x,y,o)  find  i|»(x,y,z),  where  the  field  is  propagating 
in  the  +z  direction.  Thus  for  this  problem,  we  define  S to  be  the  z=0 
plane  closed  by  an  infinite  radius  surface  such  that  V is  the  z>0  half 
sphere. 

In  (105),  G is  a suitable  Green's  function  for  the  Helmholtz 
equation  or  since  the  Green's  function  is  the  response  to  a point  source 
driving  function,  we  have  for  the  Helmholtz  eqn. 

[v2  + (2u)2k2n2]  G(x,y,z,t,n,c)  85  6(x-c)6(y-n)6(z-c)  (106) 

Using  (78)  and  (106)  in  (105)  gives 


JJJ UU.n.c)  t-(2i 


u)2k2n2  G(x,y,z,C,n,t)  + 6(x-c)6(y-n)6(z-c)} 


-G(x,y,z,c,n,t)  {-(2ir)2k2n2iU,n,t))]dV 


SJJ$  [*(C,n,o)  [- 


3G(x,y»z»S vHi 


-G(x,y,z,c,n,o)  [-  a*-k»a’*}3  dCdn  (107) 


where 


ff  CP*" 

*^c*o  plane  + 
or  S' 


Figure  4.  Surface  for  Green's  function  Integral. 

♦see  justification  on  pp  379-380  In  Born  & Wolf, 

Principles  of  Optics  [16].  The  .left  side  of  (107)  simplifies  to  give 

fffy  i(t.n.c)  «(x-C,  y-ni  z-c)  dV  * ♦(x,y,z)  (108) 

so  that  (107)  becomes 

0>(x,y,z)  = -ffs,  (i  ||  - 6 dedn  009 

Or  if  the  Green's  function  is  chosen  so  that 


(109) 


G(x,y,z,£,n,o)  = 0 


then  (109)  becomes 


;<*,*,*)  . -ffiu. n.o) 


(110) 


(111) 


(111)  is  the  desired  result.  It  gives  the  field  at  Z in  terms  of  the 
field  at  Z=0  and  the  Green's  function  G(x,y,z,t.n»c) , where  G must 
satisfy  (106)  and  (110).  Thus  it  remains  to  find  the  Green's  function. 

Fortunately,  the  Green's  function  for  the  Helmholtz  equation  with 
the  given  boundary  conditions  above  has  been  found  already  and  is  given  by  [26] 

j2*kn|f-r'|  ej27rkn|f-r| 

G(x,y ,z,€,n»0  * ~!W|ETT  + RF-FT  * (112) 


- * - i - 


where  |r-r'|  = + (y-n)2  + (z-c)‘ 

and  |r-r"  | = J(x-s)2  + (y-n)2  + (z+c)‘ 


Clearly  (112)  satisfies  the  boundary  condition  in  (110)  and  can  be 
shown  to  satisfy  the  Inhomogeneous  Helmholtz  equation 

[v2  + 4*2k2n2]G  = 6(x-U  S(yrn)  «(z-c)  - 5(x-£)  6(y-n)  6(z+c) 

(113) 

This  result  (113)  which  is  a modification  to  (106)  does  not  effect 
the  results  in  (107)  and  (108)  because  the  second  delta  function  is  out- 
side the  volume  of  Integration  V-,  and  thus  contributes  nothing  to  the 
Integral. 


Further,  it  is  possible  to  show  that  from  (112) 


3G(x,v,z.c.tuc)  . JL.  eJ2*knC(x-^+  (y-n)2  + Z2]^ 
t=o  3Z  2w[(x-02  + (y-n)2  + z2]^ 

(114) 

Then  using  (114)  in  (111)  we  have  the  solution  1>(x,y,z) 

(f  - * ( J2*kn[(x-s)2  + (y-n)2  + z2J** 


i(x,y,z)  = ^t.n.o)-^! 


2i.[(x-e)z  + (y-n)2  + zz] 


where 


i(x.y.z)  = jjdedn  g(x-t,y-n,z)  *(t,n,o) 


- < j2«kn[x2+y2+z2] 

g(x,y,z)  = --§-  ? 

32  2v[x2+y2+z2]K* 


(115) 


(116) 


(117) 


Equations  (116)  and  (117)  are  the  desired  result  and  def.ine  the  Ray lei gh- 
Sommerfeld  Integral.  If  we  perform  the  derivative  indicated  in  (117) 


rr 


2 A 


where  s ■ [x*  + y£  + 2*]  then 


a 

~n 


J2*kns  ^ ^ jknz  e 


,j2*kns 


2ws 


For  optical  frequencies  ^ Is  of  the  order  10  ®. 


ranges  less  that  10®  meters,  the  last  term  In  brackets  In  (118)  Is  much 
less  than  unity  and  is  neglected.  Further  for  s » X,  ks  » 1 and  the 
second  term  in  brackets  Is  negligible.  Thus  we  have  the  modified 
Rayleigh-Sommerfeld  result. 


i(x,y,z)  ■ 

/d"ln  -jr. 

f*J2’kSrl 

cos(z.r)  i(e,n,o) 
: — r 

isotropic,  weakly  inhomogeneous,  weakly  ab- 
sorbing medium  (r  » x) 


(119) 


1 2 2 — r 

where:  cos(z,r)  = y = cosine  between  z and  r and  r = V(x-E)  + (y-n)  +z 


Equation  (119)  is  further  modified  by  the  Fresnel  approximation 
which  holds  for  source  and  observation  points  close  to  the  z axis,  i.e. 


Then  the  following  approximations  fol  low 

cos  (z,r)  = 1 

r = d for  amplitude  term 
2 2 

r = d + ■ + (yr.nl...  + ...  for  phase  term 

Using  (121)  in  (119)  gives  the  Fresnel  Approximation 


(121) 


i>(x,y,d)  = j 

d£dj 

/ . j2nknd\ 

1.3  ) 

j2irkn  [(x-c)2  + (y-n)2] 
e d 

• i|>U.n.o) 

(122) 


Examination  of  n,  which  is  a slowly  varying  function  of  spatial 
coordinates,  yields  the  following 

n = c-^yl  = c^u(e  + |^~)  = ^|c2[u(e+  ^)]  (123) 

If  assume  that  the  magnetic  permeability  y and  the  conductivity 
o are  essentially  constant  (not  dependent  on  spatial  coordinates),  then 
the  only  dependence  on  spatial  coordinates  is  in  e.  Then  we  can  define 

y = v»o  — permeability  of  free  space 
a - constant 

e = c0  + el 

eq  = free-space  permittivity 

e,  = spatially  dependent  permittivity  due  to  inhomogeneous 
atmosphere 

Further,  for  a weakly  Inhomogeneous  medium 
C1  K<  co 


Now 


n ■ [v0(e0+ej  + ja/2nv)] 

■ V‘2  vo('  ♦ r * i 5^7> 

0 o 


(124) 
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T 


or  since  c2  • — - — (124) 


reduces  to 


Bt  + I^+J^) 


(125) 


which  is  of  the  form  1 ♦ C,  where  C is  a complex  number  with  magni- 
tude much  less  than  one.  The  square  root  is  approximated  by 


\1  + C * 1 + C 
Using  (126)  in  (125)  yields 


(126) 


n a 1 + — + j 


cot*ric  t 
Arftttertt/r 


•wveo  AW  OCX 


(127) 


If  we  substitute  (127)  into  the  Modified  Rayleigh-Sommerfeld  integral 
(119)  according  to  the  following: 

A 

-for  the  multiplicative  amplitude  term  in  (119)  use  n -x.  nQ  ■ 1 

A 

-for  the  phase  term  involving  n^  use  (127) 

Then  • • n • 

ft  / j2*kr\  ^ j2"kclr 
♦(x.y.z)  * JJdcdn  \-jxf-)  e ve°  e e0  cos(z.r)  4»(e,n,o) 


(128) 


or  if  we  define 


o = = atmospheric  extinction  coefficient 


2irk£^r 


e(r)  = = turbulence  induced  phase  shift  (recall  e,  is  also 

Go  a function  of  r) 


;<x,»,2)  - jfdtdn  e'  * tJ9<r)  cos(Z.r)  iU.n.o) 

a v jxr  • I 029) 

Equation  (129)  is  the  modified  Rayleigh-Sommerfeld  integral  extended  to 


include  atmospheric  turbulence  and  absorption.  -It  gives'  the  complex 


IV.  IMAGING  THROUGH  THE  ATMOSPHERE 


In  this  section,  the  results  of  Section  3.0  are  applied  to  the 
Imaging  sensor  problem. 

Other  Investigators  have  looked  at  various  portions  of  this  problem. 
Lutomirskl  and  Yura  [13]  give  the  development  for  the  extended  Huygen- 
Fresnel  formula.  In  [27]  Yura  develops  a general  expression  for  the 
mutual  coherence  function  of  a finite  optical  beam  propagating  In  a 
turbulent  medium,  and  in  [28]  extends  the  results  to  a sea  water  medium 
by  including  an  absorption  term. 

A later  paper  by  Lutomirski  and  Yura  [12]  goes  even  further  In 
bringing  together  all  the  necessary  ingredients  for  a comprehensive 
theory  on  the  effects  of  the  atmosphere  on  imaging  of  real  objects.  The 
development  is,  however,  still  monochromatic;  and  even  though  an 
expression  is  developed  for  modulation  contrast,  Its  usefulness  is 
limited  by  choosing  a sinusoidal  object  function  as  an  example.  Thus 
the  contrast  found  is  really  MTF  (i.e.  contrast  transfer  at  one  spatial 
frequency). 

A similar  development  for  the  MTF  through  sea  water  is  given  by 
Ishimaru  [11]. 

As  with  most  journal  papers,  the  above  references  omit  considerable 
detail  in  their  technical  development.  Thus  one  goal  of  the  present 
effort  was  to  present  a comprehensive  theory  with  essentially  all  the 
steps  included. 
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The  extended  Huygen-Fresnel  Integral,  eqn.  (68),  gives  the  complex 
scalar  optical  amplitude  at  z in  terms  of  the  field  at  z = 0 and  the 
atmospheric  point  response  function.  Next,  irradiance  levels  are  com- 
puted from  the  magnitude-squared  of  the  complex  amplitude  function.  From 
this  can  be  found  the  point-spread  function  and  modulation  transfer 
function  (MTF).  Then  by  including  the  effects  of  ambient  irradiance  it 
is  possible  to  compute  contrast  levels  at  both  the  object  and  observation 
planes.  Further,  since  the  approach  for  computing  contrast  has  been 
developed  in  terms  of  spatial  amplitude  functions,  the  dependence  on 
target  s!ze  and  range  is  included  in  the  result. 

In  (130)  it  is  pointed  out  that  e(r)  is  a random  phase  perturbation 
caused  by  inhomogeneities  in  the  atmospheric  refractive  index  along  the 
path  between  source  and  observation  points.  As  a random  phenomenon, 
quantitative  performance  measures  based  on  (130)  must  be  ensemble 
averaged. 

First  order  statistics  of  e(r)  - recall  from  the  definition  of  e(r) 


e(r)  = — — 1 r = e(x,y,d,£,n) 
eo 


where  is  the  random  quantity.  Further  is  a small  perturbation  about 
unity  and  may  be  assumed  to  be  zero  mean.  Thus  the  expected  value 


<e(r)>  = 0 


(131) 


Second  order  statistics  of  ifr(x,y,d)  - since  all  photodetectors 
respond  to  radiance  averaged  over  many  wavelengths,  our  interest  is  in  the 


quantity  (for  long  viewing  times) 

E(x,y,d)  = 0(x,y,d)  i|<*(x,y,d)  > 


(132) 


Then  from  (130)  for  the  Fresnel  approximation 

E(x,y,d)  » jfff  dcdndc'dn'  +(C.n,o)  ♦♦(c'.n'.o) 

j§^{(x-C)2  ♦ (y-n)2  - (x-c* )2  - (y  -ri,)2)/  j[e(r)  - e(r*)K 
e \e  / 


(133) 


where:  r s ^(x-T)2  ♦ (y  -n)2  + d2 

r'  ■ ^(x  - t')2  + (y  V)2  + d2 

For  an  incoherent  source 

♦ U.n.o)  **U'.n',o)  ■ | ^(C »n.o) | 2 6(C-C',  n-n') 
Using  (134)  in  (133)  gives 


-ad 

E(x,y,d)  = ? 

\<r  > 

j d£dn  E(s,n.o) 

J Irradiance  transfer  through  the  atmosphere 

(134) 


(135) 


where  E(c,n,o)  = |i|»U,n,o)n 

Equation  (135)  is  essentially  Huygen's  Principle  with  application  to 
Irradiance.  That  is,  for  the  unfocused  incoherent  source  at  U.n.o), 
the  irradiance  at  (x,y,d)  is  simply  a summation  of  spherical  wavelets 
emanating  from  all  source  points  at  approximate  distance  d and  attenuated 
by  the  atmosphere.  The  average  Irradiance  of  such  a wave  is  unaffected 
by  atmospheric  turbulence.  This  result  is  particularly  interesting  in 
light  of  the  fact  that  turbulence  very  definitely  affects  a focused  wave- 
front.  Further,  it  will  affect  the  irradiance  from  a source  that  exhibits 
partial  coherence  spatially. 

Note  that  the  ensemble  average  in  (132)  implies  a viewing  time  long 
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enough  that  the  observed  irradiance  Is  accurately  represented  by  an  ensemble 


average.  Further,  the  interchanging  of  time  and  ensemble  averages  implies 
spatial  ergodicity  for  the  fluctuating  refractive  index  term  which  further 
says  that  the  atmospheric  refractive  index  is  a shift-invariant  random 
quantity.  For  the  Fresnel  approximation,  the  assumption  of  shift 
Invariance  Implies  that  the  total  object-observation  volume  consisting  of 
a long  narrow  cylindrical  region  between  z=0,  and  z = d is  a region  of 
isoplanatism  ( i . e. , shift-invariant). 

Most  references  which  develop  equations  showing  the  effect  of  the 
atmosphere  on  optical  signals  assume  that  the  object-observation  volume 
is  essentially  isoplanatic  {[12],  [27]  - [29]  and  others.} 

A.  Simple  Lens  - Atmosphere  Imaging  System 

Next  we  consider  the  effects  of  the  atmosphere  on  optical  signals 
received  by  a simple  imaging  system.  Consider  the  system  shown  in  Figure 


Figure  6.  Simple  atmosphere-lens  Imaging  system. 


It  is  desired  to  find  ^.v,^),  the  optical  complex  amplitude  in  the  lens 

A 

image  plane,  in  terms  of  the  input  amplitude  <l<U,n,o)  and  parameters 
describing  the  lens  and  atmospheric  path. 


Starting  with  the  extended*  modified  Rayleigh-Sommerfeld  integral 
given  by  (129)  we  have 

-2J  J9(r)  /J2.krs  . 


*(x.y.n) 


dcdn  e 


Kr) 


♦U.n.o)  (136) 


where  we  have  used  cos(z,r)  « 1 for  the  lens  diameter,  DQ,  much  smaller 
than  d. 

A 

If  the  lens  transmission  function  is  represented  by  t(x,y)  then  we 

A A 

have  at  the  lens  output  plane  an  optical  signal  given  by  t(x,y)  ii(x,y,d). 
Further,  we  assume  this  optical  signal  propagates  to  the  image  plane  under 
vacuum  conditions  (i.e.,  the  path  between  lens  and  image  plane  Is  linear, 
isotropic,  non-absorbing,  source-free  and  homogeneous).  Thus  propagation 
from  the  lens  to  image  plane  Is  governed  by  the  vacuum,  modified  Rayleigh- 
Soronerfeld  integral. 

/ j2ukp\ 

*(u,v,f')  * jjdjrdy  (i^p- — J cos(z,p)t(x,y)<J»(x,y,d)  (137) 
(Image  plane  complex  amplitude] 

where:  p = ^ (u  - x)^  + (v  - y)2  + (f')^ 

and  for  f*  » ^ (u  - x)2  + (v  - y)2  ; cos(z,p)  = 1 

further,  as  shown  in  Appendix  D,  for  a thin  lens 

-j*k(x2  + y2) 


t(x,y)  = p(x,y)  e 


where  p(x,y)  is  the  pupil  function  given  by 

j2w(n,-l )d  k -a 

P(x.y)  * w(x,y)  e . w(*,.4 )C  ? 2 Do  2 

; x + y£  < (~j) 

2 2 °n  2 

* 0 i XZ  + / > (-y) 


(138) 


where:  nL  is  the  lens  refractive  index,  is  a constant,  6Q  is  the 
lens  thickness  on  its  optical  axis  and  t(x,y)  gives  the  additional 
phase  shift  through  the  lens  plus  includes  windowing  effects  from  the 
lens  aperture  function,  w(x,y).  Using  (136)  and  (139)  in  (137)  gives 


Equation  (140)  is  the  extended,  modified  Rayleigh-Sommerfeld  integral 

which  gives  the  image  plane  complex  optical  amplitude  ^(u,v,f‘ ) in  terms 

of  the  object  plane  amplitude  »|>U,n,o)  as  viewed  through  a slightly 

inhomogeneous,  slightly  absorbing  atmospheric  path. 

Next  we  consider  (140)  under  the  Fresnel  approximation,  i.e.,  the 

following  apply 

for  amplitude  terms 

1 

7 = \fu-02  + (y-n)2  + d2  - i 

1 

} - VTu-x)2  + (v-y)2  + (f1)2  - },  « ! 

ar  ad 

2 — 2 

e * e absorption  term 

for  phase  terms 

r-  ljU-tf  * (y-n)2  ♦ d2  . d ♦ ^ 

p - * <«-y)2  ♦ <n2  ■ f * ^ ^ 


(141) 


(142) 
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Then  (140)  Is  approximated  by 


j2«kd  iZ.kfff  (( 

*(u,v,f')  « e £ e e JJ  dxdyjydcdn 

je(x,y,d,e,n)  jS-jCU-O2  + (y-n)23  . 3fr{(u-x)*  + (v-y)‘3 

0 ♦U.n.o)  e 


■p(x,y)  e 


_e 

jxd 

-Jxk(x2+y2) 


iwkr 

)f, 

w 


043) 


To  shorten  the  notation  we  make  use  of  the  vector  two-dimensional  forms 

x * (x,y) 

?=U,n)  044) 

U = (u,v) 

Also,  making  use  of  (139)  and  the  constant  phase  terms  in  (143)  we  re- 
write (143)  as 


-ad/2  jA 
<l>(u,f' ) = e e 


// di  //dE 


e 

j4|S-I|3  j*|0-x|2 

e e w(x)  e 


je(x.l.d)  A 

-j 


*(*,o) 

nkx-X 


jXd 


jxf' 


(145) 


where:  A = 2itk[(d+f')  + (n^-l)dQ] 

and  W(x)  = 1 ; |x|  < 0Q/2 

0 ; |xl  > D0/2 


Equation  (145)  gives  the  Fresnel  approximation  for  the  image  plane  optical 

A A 

field  *(u,f)  in  terms  of  an  object  field  ij>(t.o)  imaged  through  the 
atmosphere,  including  absorption.  Recall  that  for  the  model  as  developed 
to  this  point  e(x,c,d)  is  the  random  phase  induced  by  small  fluctuations 
in  refractive  index  along  the  optical  path,  and  is  the  only  random 


57 


quantity  in  (145).  Strictly  speaking  the  atmospheric  extinction 
coefficient,  a,  is  also  a slowly  varying  random  function  of  the  optical 
path. 

B.  Point  Spread  Function  and  Modulation  Transfer  Function 
To  find  the  irradiance  in  the  image  plane  we  must  find 


E(u,f)  = <*(u,f)  **(u,f')  > 
Using  (145)  in  (146) 


(146) 


E(C,f)  = |~2-2  l/di  ||di,j|  d?  jjdc-  i(e,o)  i*(r,o) 


xV 


A|X-c|2-|X'-5'|2}  j^ri  |u-x|2-|u-^'  |21 
*W(x)W(x' ) e d ef 


-j%x-j?-x'*x')  j[e(x,c,d)  - e(f  .t'.d)] 

•e  e 


(147) 


J 


Again  for  an  incoherent  source  we  have 


*(5.o)  **(?’. o)  = |i|<(e,o)r  6(S-C') 


(148) 


which  physically  says  that  the  irradiance  at  a given  point  depends  only 
on  the  complex  amplitude  at  that  point.  The  result  of  (148)  on  (147)  is 
that  only  contributes  to  the  integral  and  (147)  reduces  to 

-ad 


E(u,f ' ) =■  4"d'^f,2  E^»°)  «(*)  "(*') 

J2a<|x-*l2-li,-cl2>  jf^|u-x|2-|u-x'|2}  -j^f-x-x'*?')  (149J 

• C 6 C 

j[e(r.c,d)  - eW.c.d)] 


< 


> 


Or  we  can  write 


E(u.f')  = Jjd£  E(c,tt)  G^G.f'.C.o) 


(150) 
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A 


where 


^(u.f.c.o)  mJr^  ff **([**' 

. j[e(x.l.d )-•(*', ?.d)].- 
w(x)w(5‘)  / 


JIJ{l?-«l2-|a,-C|2>  ji!k|u-x|2-|u-x'|2>  -jtyx-x-f'-i1) 

• P p p 

u»<  ('*'} 

6^  is  the  /v  atmospheric- lens  system  Green's  function.  The  expression  in 

(151)  can  be  simplified  by  recognizing  the  effect  of  certain  functions  in 

(151)  and  by  combining  the  exponential  terms,  excluding  those  in  the 

ensemble  brackets  < • > • 

We  have  the  combined  exponents 

j.k /JMl.2 . (isz) 


source-observation  volume  (i.e.,  shift  invariant),  the  MCF  does  not 
depend  on  tne  source  coordinates  \ and  further  depends  only  on  the 
difference  vector  (x-x').  Thus  we  have  for  (150)  and  (151) 


E(u,f)  = [fde  E(e,o)  s(e  + f u) 

Output  irradiance  as  of  PSF  and 

input  irradiance 

(154) 


where  S(e)  is  the  atmosphere- lens  system  point  spread  function  (PSF) 
and  is  given  by: 


-ad  <■ 

s («) 

xVf,z  J 

Jdijjdx'  w(x)w(x‘)  MCF(x-i') 

_j2ffk(x-x').e 

•e* 

Point  Spread  Function  For  Atmosphere-Lens  System 

(155) 


Further,  if  we  make  the  following  change  of  variable 


p * t-r 


(156) 


then  X a R + ; X'  = R - jj- 

Substituting  (156)  into  (155)  yields 


! 

jdp  JJdR  w(R  + |)  w(R  - |) 

•MCF(p)  e d 

Alternate  Form  of  Point  Spread  Function 

(157) 


Now  eqns.  (154)  and  (157)  give  the  relationship  between  source  irradiance 
and  image  irradiance  as  functions  of  spatial  coordinates*.  If  we  examine 
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1 


-j2*K*u 


the  relationship  in  the  spatial  frequency  domain,  then  from 

- . (f  -j2*iC-u 

E(K,f)=jjdO  E(u.f)  e 

and  using  (154)  in  (158) 

E(K.f*)  = jjdu  ffdc  Etf.o)  S[t  + f 8) 

or  with  the  following  change  of  variable 

V • l * f u 

u = J-(C'-t) 
du  = f’/d  d«* 

(159)  becomes 
E(K, 


(158) 


(159) 


i.n  -fjjai 


.1  - . j 

• S(c')  e d JJ d*'  EU*o)  e 


2irf  ’K-l 


Then  from  the  Fourier  transform  definition  (160)  is  given  by 


£(ic,f)  - £ s(*r)E  (-*£-.  c) 


(160) 


(161) 


Spatial  Frequency  Domain 


where:  K is  a two  dimensional  spatial  frequency  vector  (cycles/meter) 

A — 

S(K)  is  the  Fourier  transform  of  the  point  spread  function  (eqn. 157) 
and  is  by  definition  the  total  unnormalized  atmosphere-lens  system  optical 
transfer  function  (OTF). 

The  interpretation  of  equation  (161)  is  as  follows.  If  the  object 

A 

irradiancc  function  is  known(i.e.  E(£,o))  then  the  image  Irradiance  function 

A A » 

E(x.f')  is  related  to  E(e,o)  and  the  atmosphere-lens  transfer  function  by 
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equation  (161)  In  terms  of  the  spatial  frequency  content  of  the  ohject  and 
Image.  The  negative  argument  in  the  last  term  of  (161)  Is  a result  of 
image  inversion  by  the  lens.  Under  the  assumptions  of  circular  symmetry 
for  the  lens  and  isoplanicity  for  the  atmospheric  path,  the  only  effect 
of  the  negative  sign  is  a reversal  in  reference  direction  for  the  image. 
The  spatial  frequency  functions  will  be  unaffected. 


Examination  of  the  inner  integral  in  (157)  yields 


Q(5)  = jj dR  w(R  + |)  w(R  - |) 


(162) 


where  the  integral  is  over  the  lens  aperture  and  (163)  is  shown  to  be 
[32] 

4' [»-«)■ 

=0  J |p|  > D0  (163) 

. 2 
The  value  of  Q at  |p|  = 0 is  nDQ  /4 

so  that  by  normalizing  Q we  have  the  lens  modulation  transfer  function 


(164) 


and  equation  (157)  can  be  written  as 


0,d  *°„2  ff 

qi  -rJ)dS 


-j2trkp»£ 

M(p)  MCF(p)  e d 


Now  since  the  optical  transfer  function  is  the  Fourier  transorm,  S(K),  of 
S(£)  and  letting 


K > s JL 

d xd 


p - Kxd 
dp  = Xd  dK 


(166) 


Using  (166)  in  (165)  gives 


O „ 

*Dn e ((  _ _ „ -j2wK*e  , v 

S(c)  = „ JdK  M(KXd)  MCF(KXd)  e (167) 

4xJdf,z 


o uu 

. - *Dne 

OTF  = S(K)  = 5.  M(KXd)  MCF(KXd) 

4x3df,Z 


068) 


Since  the  system  MTF  is  the  normalized  magnitude  of  the  OTF  we  have 
finally  that  the  atmosphere- lens  system  MTF  is 


M(KXd)  MCF(KXd) 
MCFTo) 


069) 


Thus  determination  of  the  atmospheric  MTF  consists  of  determining  the 

mutual  coherence  function,  MCF. 

C.  Atmospheric  Wave  Structure  Function 


The  approach  given  here  is  based  or  work,  by  Fried  [15]  and  considers 
a slmplt  lens  system  illuminated  by  an  atmospherically  degraded  plane 
wavefront.  The  result  gives  MTF  for  the  system  in  terms  of  the  atmospheric 
wave  structure  function. 

Consider  a simple  imaging  system  as  shown  in  Figure  7.  Then  for  a 
wave  incident  on  the  lens  which  deviates  little  in  phase  and  amplitude 

a 

from  a plane  wave  and  designated  as  <|>0(x,y,d)  the  resulting  wave  in  the 

focal  plane  of  the  lens  is  shown  to  be  [16]  approximated  by 

* C ( -Jp*  ( xx ' + yy ' ) 

*(u,v,z)  = A JJdxdy  *0(x ,y,d)e  x (170) 

where  X Is  the  optical  wavelength 
f is  the  lens  focal  length 
and  k * i is  the  wavenumber 


*r**i 


X 


Figure  7.  Focal  plane  field  for  a simple  lens 

If  we  use  the  shortened  notation  p'  = (u,v)  and  p * (x,y)  then  (170) 

Is 

- ( ( - . p 1 ) 

♦(?'  ,2)  = A JJ  dp'  <i0(p*  ,0)e  (171) 

The  MTF  of  the  total  system  (optics  plus  atmosphere)  is  simply  the 
normalized  two-dimensional  spatial  Fourier  transform  of  the  intensity  of 
the  image.  That  is 

MTF(K)  = B jj  dp  *(r  )J*(p*  (172) 

where  K = (e  ,B  ) spatial  frequency  in  x and  y directions 
a y 

B is  a normalizing  constant  which  makes  MTF(O)  * 1. 

Next,  expressing  the  input  wave  as 

i(x,y,0)  ■ W(x,y,d)e^w,y)  + (173) 

where  W(p)  = H(x,y,d)  * 1;  |p|  <.  0o/2 

= 0;  otherwise 
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is  the  lens  pupil  function.  i(p)  and  $(p)  are  the  log-amplitude  and 
phase  of  the  incident  wave  respectively,  i and  $ will  be  random  variables 
because  of  atmospheric  refractive  index  fluctuations;  therefore,  MTF(K) 
will  be  a random  variable  also.  Thus  we  must  describe  the  MTF  in  terms 
of  statistical  averages. 

For  the  long-time  exposure  case  Fried  shows  the  average  MTF  to  be 


where 


<MTF(K)>  LE  = MTF0(K) 


- \ D(xfK) 
e 


(174) 


MTFq(K)  is  the  lens  MTF 

O(xfK)  is  the  wave  structure  function  for  the  atmospheric 
path  and  in  terms  of  (173) 


MTFq(K)  = 


A2B 


If 


dp  W(p  - AfK)W(p) 


D(XfK)  = [♦(p)  - + (p  - XfK)]2 

+ U(p)  - *(p  - XfK)]2. 

For  this  case  the  effect  of  the  atmosphere  is  separable  from  that  of  the 
lens  and  is  given  by: 


<MTF(K‘)>  ATM0SPHERE  = 


LONG  TIME  EXPOSURE 


(175) 


where  K*  has  units  of  cycles  per  radian  field-of-view;  K'  = fK.  Thus 
(175)  gives  the  MTF  of  the  atmosphere  in  terms  which  are  independent  of 
the  lens  parameters. 


65 


- r 


For  the  short-exposure  case  the  result  is 

<MT m>  SE  . ■*,«-«.-  i 6,1,6  i ’ - Hf ) ’1 

076) 

where  "b"  is  a constant  between  for  the  far-field  and  unity  for  the 
near  field.  In  this  case  the  effect  due  to  the  atmosphere  cannot  be 
expressed  independently  of  the  lens  parameters,  i.e. 


<MTF(R')^TM0SpHERE  = exp 

- D(xK' ) j 

1 

" b 

/ - \ 1/3,il 
© ] 

SHORT  EXPOSURE  ' 

In  a later  paper,  Metheny  and  Philbrick  [30]  present  an  argument 
based  on  the  Taylor  "frozen  flow"  hypothesis  for  separating  the  optical 
transfer  function  for  exposure  times  considerably  shorter  than  those 
required  for  an  ensemble  average  as  was  needed  in  (175). 

Further  insight  into  the  results  obtained  in  this  section  is  gained 

A 

from  the  following  considerations.  The  wave  structure  function  D(r)  has 
been  shown  theoretically  and  experimentally  [1]  and  [31]  to  have  the 
relationship 

D(r)  = Ar5/3  * 6.88  5/3  (178) 

where:  r^  = (6.88/A)3^5  and  the  factor  6.88  is  arbitrarily  chosen  to 
make  the  knee  of  the  resolution  function  occur  at  Tq  * Dc(see  Fried  [15] 
for  a discussion).  The  factor  "A"  is  a constant  which  depends  on  the 
propagation  path  length,  the  wavelength,  the  strength  of.  the  turbulence 
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and  the  nature  of  the  unperturbed  wavefront.  Equation  (178)  has  been 
shown  to  be  valid  for  both  plane  and  spherical  wavefronts.  Substituting 
(178)  into  (175)  and  (177)  we  have  finally  the  MTF  due  to  the  atmosphere 
in  terms  of  the  parameter  rg. 


-3.44 


<MTF(K')>  LE 

ATMOSPHERE 


(?) 


(179) 


* aa  /*K'\  5/3  ( , 1 /XK'\  1/3) 

-3-44(^)  p-F  (nr)  ) 


<MTF(K')>  se 

ATMOSPHERE 


(180) 


D.  Contrast  Based  on  Em  Theoi 


The  total  irradiance  in  the  image  plane  u = (u,v)  is  the  received 
object  irradiance  E given  by  equation  (154)  plus  ambient  irradiance  E . 


ET0T  = E + Ea 


(181) 


From  which  we  can  define  universal  contrast  at  the  image  plane 


A {ET0t} 


MAX  " TTOTIMIN 


{ etot}k 


A EMAX  ' EMIN 


emin  + eamin 


Then  using  (154)  in  (183)  gives 


c „ [ l\dj  ig,o)  S(j  ♦ fru)]^  -'[{/  dC  E(6.o)  S(l  ♦ fo)], 

[/|dc  E ( C »o ) S[j  + ^Tti)]  MIN  * EjMIN 


IMAGE  CONTRAST  FOR  ATMOSPHERE-LENS  SYSTEM 


Equation  (184)  gives  the  contrast  In  the  image  of  an  atmosphere- lens 
Imaging  system  in  terms  of  the  atmosphere-lens  point-spread  function 
S(l),  the  ambient  irradiance  Ea  and  the  object  Irradiance  function  E(£,o). 
Inherent  In  the  expression  are  dependence  on  atmospheric  extinction  and  on 

A 

object  dimensions.  Further  analysis  requires  that  E(£,o)  be  specified  and 
that  S(t)  be  known. 

An  alternate  expression  for  ^ is  obtained  In  terms  of  spatial 
frequency  functions 


r a _ j2*K-t.  - 

{J_.  E(K,f')e  dK>  - { f_m  E(K,f' )e 


j2*K*S_ 
dK) 


*1  *R 


MIN 


- ~ — MT. 

{ J E(K,f ' )e  dK) 


(185) 


MIN 


where  E(K,f')  is  the  image  spatial  transform  as  given  by  equations  (161) 
and  (168). 

If  E(£,o)  is  a one-dimensional  periodic  function  such  as  would  be 
represented  by  a bar  chart  (see  Figure  below)  then 


Figure  8.  Bar  pattern. 


E(*,o)  = £ *r 


n=-® 


S(K<  • v 


(186) 
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where  X0  is  the  object  period  in  the  X-direction.  The  object  Is  then 

represented  by  its  Fourier  series 
« j2irnX 

E(C.o)  = ^ En  e *o  (187) 

n=-« 


and  in  the  image  plane  from  (187)  and  (161),  for  sinusoids 

s2irf 'n< 

'o  (188) 

n=-« 

Now  if  the  contrast  of  the  bars  is  computed  using  X = o as  the  maxi- 


•  ,2irf'n 

-E  f K ^ 


mum  irradiance  coordinate  and  x = dxo/2f'  as  the  minimum  irradiance 
coordinate  in  the  image  plane  then 

n=-“  ° 

and 

^•f'U‘E  f'-""  5(o) 


(189) 


n=-« 


Then 


'1  »R 


?E  K 


n=-®> 


n odd  only 


E *(£K*  f* 


n=-» 


(190) 


UNIVERSAL  CONTRAST  PERIODIC  BAR  CHART 
THROUGH  ATMOSPHERE-LENS 


Equation  (190)  is  the  desired  result  and  gives  contrast  in  terms  of 

A 

the  atmosphere-lens  optical  transfer  function  S the  ambient  irradiance 

* A 

Efl  and  the  Fourier  coefficients,  En,  for  the  bar  pattern.  The  expression 


y~ 
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also  includes  atmospheric  absorption. 


In  terms  of  the  modulation  contrast,  R,  we  have 


j»R 


EMAX  “ EMIN 


EMAX  + EMIN 


Thus  for  the  bar  chart  image 


V.  TEST  AND  MEASUREMENT  METHODS 

A.  Visibility  and  Contraat 

For  the  contrast  reduction  expression 


*» 


V 

M ' ® 

nr 


(193) 


complete  specification  requires  measurement  of  the  background  radiances 
Nfl'  and  Nr'  at  ranges  of  0 and  R.  meters  in  addition  to  "visibility" 
measurements.  Since  the  visibility  range  Rv  is  defined  to  be  the  range 
at  which  contrast  is  reduced  to  2 percent  of  CQ,  then  Ry  lends  itself 
well  to  measurement  with  a telephotometer.  Using  measured  values  of 
Ry,  Nq'  and  NR' , we  can  compute  the  value  of  the  extinction  coefficient, 
i.e. 


(194) 


Once  Oq  is  known,  contrast  measurements  may  be  taken  at  various 
ranges  and  compared  with  the  predicted  values  of  (193).  The  comparison 
is  in  fact  a comparison  of  R the  optical  slant  range  with  R,  the  true 
slant  range. 

If  the  measurement  situation  is  for  a horizontal  path  with  an  object 
being  measured  against  a sky  background  then  NR'  = N ’ and  it  is  usually 
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assumed  that  R = R so  that  measurement  of  contrast  versus  range  will 
give  all  the  pertinent  information.  The  result  should  be  as  shown  in 
Figure  9.  Extinction  coefficient  is  then  easily  obtained  from  the 
slope  of  the  curve  while  visibility  range  is  that  range  at  which 

Vco  = -02- 

For  the  case  of  (193)  where  the  path  is  not  horizontal  against  a 
sky  background,  the  effects  of  N^'  and  R must  be  considered.  Since 


Figure  9.  Visibility  and  extinction  coefficients  from 
contrast  measurements. 


Nr'  is  the  background  radiance  at  range  R and  since  all  objects  (bright 
or  dark)  tend  to  approach  the  air  light  radiance  at  long  ranges,  then 
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NR'  will  decrease  with  R for  a bright  background  and  Increase  with  R 
for  a dark  background.  Thus  an  object  against  a background  can  have  its 
contrast  reduced  either  more  rapidly  or  less  rapidly  for  the  sky-to- 
ground  case  than  for  the  horizontal  case.  With  respect  to  Figure  9, 
this  effect  would  be  realized  as  a change  In  slope  of  the  curve.  In 
fact.  It  is  not  unreasonable  to  expect  that  small  nonlinearities  may 
exist  in  the  NR'  dependence  on  R.  The  optical  slant  range  R is  expected 
to  be  a gradually  decreasing  function  of  R since  it  is  a measure  of  the 
optical  path  which  depends  on  refractive  index.  The  effect  of  Figure  9 
is  to  rotate  the  curve  slightly  .counter-clockwise  about  Cr/Cq  = 1-0. 

A number  of  precautions  and  effects  must  be  considered  in  performing 
the  indicated  measurements.  First,  since  any  measuring  instrument  will 
be  dependent  on  spatial  frequency  content  of  a test  target,  this 
dependency  must  be  known  and  compensated  in  interpreting  the  measurement 
data.  This  problem,  of  course,  arises  because  of  the  size  with  range 
variation  of  what  the  measurement  instrument  sees.  Otherwise  some  zoom 
mechanism  must  be  used  in  conjunction  with  the  instrument  to  maintain 
constant  relative  image  size.  If  such  an  arrangement  is  used,  however, 
it  is  still  necessary  to  consider  the  MTF  of  the  zoom  optics  and  correct 
for  It. 

B.  Wave  Structure  Function  and  Atmospheric  Mutual  Coherence  Function 

MTF  effects  of  the  atmosphere  were  shown  In  Section  4.3  to  depend 
on  a quantity  6(r)  which  was  called  the  wave  structure  function  where  "r" 
Is  a correlation  parameter  representing  spatial  separation.  It  has  been 
shown  for  a large  range  of  values  "r"  that  the  wave  structure  function  is 
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adequately  represented  by  [33] 


D(r)  - Cn2  r2/3  (195) 

where  C 2 Is  the  refractive  index  structure  constant,  or  more  accurately, 
n 

2 

the  refractive  Index  structure  function,  since  C Is  not  a constant, 

n 

but  varies  with  time  and  space.  In  fact,  as  seen  by  curves  In  [34], 

2 

Cn  can  vary  over  several  orders  of  magnitude  during  a typical  day.  It 
Is  strongly  dependent  on  the  time  of  day  as  well  as  altitude.  The  wave 
structure  function  can  be  written  as 

D(r)  = 6A(r)  + 6$(r)  ' (196) 

A 

where:  Dft(r)  is  the  variance  of  the  amplitude  difference  fluctuations 
of  the  refractive  Index  between  two  points  spaced  "r"  apart  and  D^(r) 

A 

Is  the  similar  variance  of  the  phase  terms.  Although  DA(r)  may  be 

A 

measureJ  using  a pair  of  photodetectors  [1]  and  Ds(r)  may  be  measured 

with  an  interferometer  [35]  the  most  commonly  cited  technique  is  to 

2 

measure  a function  CT  which  is  defined  as  the  temperature  structure 

2 

function  using  temperature  profile  measurements.  From  C^.  the  refractive 

2 

index  structure  function  C is  calculated  which  is  then  related  to  the 

n 

MTF  of  the  atmosphere. 

For  optical  frequencies  a good  approximation  for  the  refractive 
index  "n"  is  given  by  [36] 

n » 1 + A(p/T)  [1  + B -|jr]  (197) 
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where  p = pressure  (millibars) 

T = temperature  (Kelvin) 
e = vapor  pressure  (millibars) 
B is  a constant,  and 


A « 77.5(10”6)[1  + 5.15(10"3)X"2  + 1 . 07 ( 1 0-4 ) X-4] 


(198) 


is  a wavelength  dependent  constant.  Higher  order  terms  in  (197)  are 
neglected.  It  can  be  shown  that  if  humidity  effects  can  be  ignored 
then 


C_2  <v,  cx2 

T2 


where  CT2r2/3  = <\[T(r1 ) - T(r2)]2^> 


099) 

(200) 


r = lri  ‘ r2| 


Examples  of  Cy  for  various  temperature  measurement  conditions 
(eg.  altitude,  time  of  day,  etc.)  are  given  in  [37].  Example  MTF 
curves  are  given  in  the  same  reference  for  various  viewing  conditions. 

The  atmosphere-lens  imaging  system  is  specified  by  either  its 
point  spread  function  (PSF)  equation  (157)  or  its  modulation  transfer 
function  (MTF),  S(^-)  in  equation  (161). 

From  (161)  we  can  solve  for  the  total  system  MTF 


d JiC^j 


(201) 


where  E(K,f')  is  the  two-dimensional  spatial  fourier  transform  of  the 
output  image  irradiance  function.  This  quantity  can  be  computed  digitally 
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from  a measurement  of  the  two-dimensional  irradiance  function  in  the 
image  plane  of  a measurement  instrument.  E^-  is  the  scaled  two- 

dimensional  spatial  fourier  transform  of  the  input  object  irradiance. 

It  can  be  measured  or  computed  from  knowledge  of  the  object  to  be  imaged 
(eg.  a bar  chart  pattern).  The  effect  of  ambient  irradiance  or  average 
irradiance  level  of  the  bar  chart  is  on  only  the  K = o term  for 
equation  (201).  Thus,  these  effects  can  be  quantized  and  accounted  for. 

The  computation  and  scaling  for  the  input  spatial  transform  is  best  done 
digitally.  Further  to  avoid  overflow  problems,  the  result  in  (193) 
should  be  computed  only  for  values  of  such  that  the  denominator  is 
larger  than  zero,  and  in  fact,  larger  than  some  system  noise  level. 

The  method  outlined  above  should  give  reasonable  accuracy  for  the 
atmosphere- lens  system  at  a given  range.  Measurements  at  different  ranges 
shorld  affect  the  MTF.  In  particular,  the  absorption  term  in  equation  (157) 
varies  exponentially  with  range  "d".  Further,  the  statistical  characteristics 
of  the  mutual  coherence  function  may  change  with  d!  The  complete  dependence 

**  rt 

on  d is  difficult  to  predict  analytically.  However,  a measurement  program 
which  includes  several  ranges  within  a total  expected  range  variation 
could  provide  data  for  experimental  prediction  of  the  total  range  effect. 

Essentially  the  experimental  program  should  consist  of  measurements 
of  image  plane  irradiance  for  computation  of  equation  (201).  The  same  bar 
chart  target  can  be  used  at  the  different  ranges;  however,  care  must  be 
exercised  so  that  the  received  image  size  Is  not  so  small  as  to  be  down 
in  the  noise  of  the  measuring  system  MTF.  An  instrument  with  adjustable 
field-of  view  should  be  capable  of  compensating  for  target  size  variations 
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over  a typical  range  variation  of  one  to  five  kilometers. 

For  determining  contrast  transfer  through  the  atmosphere  lens  system, 
equations  (184)  and  (185)  give  the  theoretical  results.  A result  more 
amenable  to  experimental  varification,  however,  is  to  use  equation  (182) 
where  the  actual  image  irradiance  levels  can  be  measured.  This  result 
can  be  compared  with  the  inherent  object  irradiance  to  get  contrast 
transfer  through  the  system. 

The  effect  of  target  size  on  contrast  transfer  is  easily  measured 
by  including  bar  charts  of  varying  size  at  each  range.  Care  must  be  taken 
to  correct  for  the  lens  effect  if  it  is  desired  to  know  the  effect  of  the 
atmospheric  path  only. 
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In  summary,  the  methods  for  measuring  effects  of  the  atmosphere  on 
an  optical  system  can  be  as  complex  as  the  theory  which  predicts  these 
effects.  A more  detailed  study  of  this  topic  should  be  done  to  fully 
investigate  the  validity  of  specific  measurement  techniques  in  particular 
as  related  to  theoretical  results.  Only  after  evaluation  of  such  an 
effort  would  it  be  practical  to  investigate  simulation  methods. 


VI.  SUMMARY  AND  CONCLUSIONS 


In  evaluating  the  usefulness  of  the  theoretical  expressions  developed 
in  this  report  in  terms  of  their  applicability  for  determining  performance 
of  imaging  seekers  it  becomes  necessary  to  examine  the  meaning  of  various 
assumptions  and  approximations  made  during  the  course  of  those  theoretical 
developments.  For  example,  some  method  for  reconciling  a monochromatic 
theory  with  a wideband,  spectrally  non-trivial  application  such  as  is  true 
for  the  visible  and  near-infrared  imaging  seeker  must  be  found.  This 
problem  has  not  been  treated  rigorously  in  the  literature,  the  contrast 
theory  as  developed  by  Duntley  and  Middleton  has  been  extended  to  give 
results  for  a finite  wavelength  band  and  was  shown  to  have  little  effect 
except  for  slant  ranges  greater  than  or  close  to  the  meteorological  range. 
Under  these  conditions  there  is  of  course  serious  doubt  that  the  sensors 
can  operate  at  all. 

If  we  nowexamine  the  implications  of  a finite  spectral  band  on  the 
results  of  Section  4.0,  thenthe  MTF  of  the  atmosphere  for  imaging  seeker 
applications  should  be  the  result.  It  is  between  equations  (75)  and  (76) 
that  a monochromatic  field  was  assumed.  The  major  advantage  of  this 
assumption  was  that  it  allowed  separation  of  the  time  and  spatial 
dependence  of  the  field.  Therefore,  a solution  $(x,y,z,v)  which  is 
not  necessarily  monochromatic  should  still  satisfy  the  scalar  Helmholtz 
wa«c  equation.  Unfortunately  the  approximate  solution  methods  take 
advent age  of  the  separability  and  in  some  cases  assume  a propagating 


plane  wave  so  that  these  methods  are  not  easily  extended  to  the  case 
of  a continuous  spectral  band  and  a general  $(x,y,z,v). 

If  a monochromatic  wave  4>(x,y,z,v0)  satisfies  the  wave  equation, 
then  a summation  of  monochromatic  waves  should  also  satisfy  the  wave 
equation,  i.e. 

♦ s Yj  ♦ , (x,y,z,vj ) (202) 

and  if  we  extend  this  concept  over  a continuous  range  of  optical 
frequencies  we  get  as  a solution  also 

♦ s / $(x,y,z,v)dv  (203) 

•'Av 

Further  development  and  evaluation  of  this  approach  are  needed.  A 
different  approach  based  on  the  theory  of  partial  coherence  is  another 
alternative  for  the  prediction  of  broadband  MTF  effects  of  the  atmosphere. 

The  development  given  in  this  section  for  including  broadband  effects 
on  the  contrast  transfer  function  is  easily  applied  to  the  imaging  sensor 
application.  The  development  given  for  broadband  effects  on  the  MTF 
due  to  atmospheric  turbulence  requires  further  work. 

In  conclusion,  further  details  of  the  experimental  program  need  to 
be  worked  out.  Additionally,  there  are  still  some  rough  spots  in  the 
theoretical  development,  which  if  solved,  could  provide  more  useful 
analytical  predictions  of  the  effects  of  range  and  target  size  on  MTF 
ano  contrast  transfer. 

The  comprehensive  rigorous  theory  which  was  developed  in  this 
effort  lays  the  groundwork  for  a more  complete  description  of  the  finite 
spectral  band  atmospheric  imaging  system. 
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APPENDIX  A:  VISIBILITY,  VISIBILITY  RANGE 
AND  METEOROLOGICAL  RANGE 


Visibility,  visibility  range  and  meteorological  range  are  three 
terms  for  the  same  quantity  which  is  defined  specifically  as  the  range 
at  which  the  apparent  contrast  of  an  object  Is  two  percent  of  its 
inherent  contrast.  Thus  "visibility"  has  units  of  distance  and  is  a 
somewhat  arbitrary  measure  of  how  "clear"  the  atmosphere  Is  on  any  given 
day  or  for  a given  set  of  meteorological  conditions.  For  the  simplest 
contrast  transfer  relationship  given  in  equation  (12)  if  Ry  Is  the 
visibility  range  then 


co  -a0^v 

= .02  = e U v 

c0 


(A-l) 


r 


1 


APPENDIX  B:  CONTRAST  DEFINITIONS 

A number  of  definitions  for  contrast  are  found  In  the  literature 
of  which  three  commonly  used  definitions  are  given  below.  Judgement  as 
to  which  Is  best  Is  left  to  the  discretion  of  the  user;  however,  the 
significant  differences  are  identified. 


'1 


A L0  ~ Lb  = ^0  . i.  .1  < r < 
L L.  '• 


lLn  " 


r A ^0  “ ^ 

c3  = TL^n-y 


r'1 

Lb 


r-1 

Lb 


; 0 £ ® 


(Universal ) 

(Television) 


5 

u 


;o  + 1 


-1  <.  C3  <_  1 (Modulation) 


Definition  has  the  property  that  for  any  given  background  radiance 
L^,  the  contrast  increases  linearly  as  the  object  radiance  increases. 

It  takes  on  values  of  -1  for  a totally  dark  object  to  +«  for  a totally 
dark  background.  Definition  C2  is  similar  to  except  that  the  abso- 
lute magnitude  sign  eliminates  any  negative  contrast  values.  The  argument 
is  that  the  important  thing  is  radiance  difference  between  object  and 
background  and  not  which  is  brighter.  Definition  eliminates  large 
numbers  and  maps  contrast  into  the  range  -1  to  +1 ; however,  the  contrast 
dependence  on  increasing  object  radiance  is  now  a non-linear  function. 

in 

All  three  contrast  definitions  are  shown  plotted^Figure  B -1  versus 
L0/Lb* 
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APPENDIX  C:  LINEAR  SYSTEMS,  POINT  SPREAD  FUNCTION, 
MODULATION  TRANSFER  FUNCTION  AND 
CONTRAST  TRANSFER  FUNCTION 


1 


Based  on  assumptions  of  a weakly  inhomogeneous,  weakly  absorbing 
atmosphere  with  enough  uniformity  that  the  target/imaging  system  can  be 
considered  to  be  in  an  isoplanatic  region  (i.e.,  the  atmospheric 
region  of  interest  can  be  assumed  to  be  statistically  shif-invariant) 
and  taking  into  account  the  size  scaling  and  reversal  of  images  from  a 
simple  lens  system,  the  atmosphere- lens  system  is  a quasi-linear,  shift- 
invariant  system  and  can  be  described  in  the  following  way. 

If  the  object  to  be  imaged  is  represented  by  its  complex  monochromatic 
spatial  amplitude  function  ii<(x,y,o),  then  the  image  spatial  amplitude 
function  is  given  by  the  convolution  r»f  g(x,y,z)  with  i)»(x,y,o). 

*(x,y,z)  = g(x,y,z)  * ^(x.y.o)  (Cl) 


i 


where  g(x,y,z)  is  the  two-dimensional  green's  function  for  the  atmosphere- 
lens  combination  and  z = d + f'  is  the  total  object-image  separation. 

See  Figure  Cl . 


r - Zr*/3  Focal  le*s6TH 
Figure  C-l.  Atmospheric-lens  imaging  system. 
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Or  schematically,  equation  Cl  can  be  represented  by  Figure  C2 


Figure  C-2.  Schematic  representative,  linear  system. 

Our  interest  is  in  the  real  image  intensity  and  not  in  the  complex 
spatial  amplitude  function.  So  we  have  the  image  irradiance  function 
of  spatial  coordinates 

E(x,y.z)  = <«|i(x,y,z)<(»*(x,y,z)>  (C2) 

where  the  brackets  denote  an  ensemble  average  and  are  necessary  because 
the  atmospheric  green's  function  is  a random  quantity. 

From  (Cl)  and  (C2) 

E(x,y,z)  = <{g(x,y,z)  * ii»(x,y,o)}*{g*(x,y,z)*<|>*(x,y,o)}  > (C3) 

where  the  superscript  * denotes  complex  conjugate.  Or 

E(x,y,z)  =<^jj|jdad0dYdn  g(x-a,  y-B,  z)  ti>(a,B,o) 

■9*(x-y,  y-n,  z)  ^Y.n.o)^  (C4) 

If  the  source  is  spatially  incoherent  then 

A A A 

,o)  V^Y.n.o)  = E(a,B,o)  6(y-ct,  n-B)  (C5) 

and  (C4)  reduces  to 


E(x,y,z)  =(}(dadg  g(x-a.y-B.o)  g*(x-a,y-B,o) 


Equation  (C6)  expresses  the  image  plane  irradiance  E(x,y,z)  in  terms 
of  the  object  plane  irradiance  E(a,B.o)  and  the  mean  square  value  of  the 
atmosphere- lens  green's  function.  Can  also  write  (C7)  as 


E(x.y.z)  * |g(x,y,z)|2  * E(x,y,o)  (C7) 

or 

E(x,y,z)  = s(x,y,z)  • E(x,y,o)  (C8) 

where:  s(x,y,z)  » |g(x,y,z)|2  is  the  atmosphere-lens  point  spread 

function.  That  is  s(x,y,z)  is  the  two  dimensional  image  irradiance 
function  at  z = z for  an  ideal  point  source  input  irradiance  function 
at  z = o. 

From  the  convolution  theorem  and  taking  the  spatial  Fourier  transform 
of  (A8)  we  get 

E(kx,  ky,  z)  = ^{s(x,y,z)}E(kx,  ky,  o)  (C9) 

where  E(kx,  ky»  o)  is  the  input  irradiance  spatial  frequency  function 
E(kx  , ky,  z)  is  the  output  image  irradiance  spation  frequency 
function 

and  ky,  k^  are  spatial  frequencies  in  x and  y direction  respectively. 

Further 

I' 

J{s(x,y,z)}  = S(kx,  ky,  z)  (CIO) 

| where  $(kx,  ky,  z)  is  the  combined  atmosphere-lens  optical  transfer 

[ function  for  monochromatic  radiation. 

The  modulation  transfer  function  for  the  combined  system  is  then 

HTFIK,.  v 2 ) - J%7k^  <cm 


i 


The  MTF  as  expressed  In  (Cl  1 ) gives  the  total  system  normalized 
response  to  sinusoidal  spatial  irradiance  as  a function  of  spatial 


frequency. 

This  result  may  be  compared  to  the  contrast  transfer  function  (CTF) 
in  which  case  CTF  is  the  normalized  total  system  response  to  a square- 
wave  spatial  irradiance  as  a function  of  spatial  frequency. 

In  terms  of  equations  (C8)  and  (C9)»  to  get  MTF  the  input  irradiance 


function  is  a pure  sinusoid,  i.e.. 


E(x,y,o)  ■ ^ (1  + cos  2n{ko>x  x + kQ>y  y}) 


(C12) 


E(kx,ky,o)  - 2°  6(kx,ky)  + 4 <5^x’lco,x’lcy”ko,y^ 


+ -f-  6(kx  + k0>x,  ky  + kQ>y) 


(C13) 


To  get  CTF,  the  input  irradiance  function  is 


E(x,y,o)  = SQ(x,y,o) 


(Cl  4) 


where  SQ(x,y,o)  is  a two  dimensional  square  wave  as  shown  in  Figure  C3. 


Figure  C-3.  General  two-dimensional  bar  pattern. 


j 
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The  two  dimensional  Fourier  spatial  transform  of  (C14)  gives 


E(kx,ky,o)  = £ Am  «(kx-nik0fX,ky-mk0>y) 


(C15) 


where  the  coefficients  A_  are  given  by 

m 


Am  = 


and  sine  6 sin 
wx 


(C16) 


(C17) 


Thus  the  square  wave  input  irradiance  function  consists  of  a summation 
of  weighted  sinusoidal  spatial  irradiance  functions.  Then  it  follows  that 
the  response  to  a square  wave  input  is  the  sum  of  the  responses  to  the 
weighted  sinusoids  comprising  the  square  wave.  Then  we  have  the  following 
relationship  between  CTF  and  MTF. 


- £»„  MTF(kx>m,kyiln.z) 


(C18) 


So  in  terms  of  the  real  image  irradiance  functions  the  following  linear 


system  concept  applies. 


E(x,y,o) 


s(x,y,z)  . 
=|g(x,y.z)|‘ 


E(x,y,z) 


The  remainder  of  this  Appendix  gives  a different  approach  to 
definition  of  MTF  and  gives  some  useful  relationships  among  several 
optical  system  performance  functions. 

Analogous  to  the  transfer  function,  H(f),  of  an  electrical  network, 
we  can  define  an  optical  transfer  function  (OTF)  which  relates  input 
and  output  spatial  frequency  content  for  an  optical  system.  The  OTF 


is  then  given  by 


I(k„,kJ 


0TF  * snyy  • A(lvVc 


je(kx,ky) 


(Cl  9) 
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where  I(k  ,k  ) = image  plane  (output)  representation  as  a function 

* y 

of  spatial  frequencies  k and  k 

x y 

0(k  ,k  ) = object  plane  (input)  representation, 
x y 

The  OTF  is  in  general  complex  so  that  A(k  ,k  ) is  its  amplitude  and 

^ y 

e(k  ,k  ) is  its  phase.  Both  are  functions  of  spatial  frequency.  There 

* y 

is  also  an  implicit  dependence  on  spectral  content  of  the  radiation. 

The  computation  and/or  measurement  of  the  optical  transfer  function 

is  accomplished  as  a time-averaged  value  since  all  practical  photosensors 

have  finite  response  times  and  are  in  fact  sensitive  to  the  optical 

Intensity,  not  the  amplitude.  However,  for  measuring  the  OTF  of  a 

slowly  time-varying  medium  such  as  the  atmosphere,  the  observation  time 

becomes  an  important  parameter.  In  other  words,  any  observation  will 

include  a large  number  of  optical  wavelengths  over  which  the  effect  Is 

averaged,  but  the  given  wavelength  train  can  be  modulated  by  the  properties 

of  the  atmosphere  at  rates  which  are  measurable. 

Since  the  spatial  phase  function  is  essentially  random  for  incoherent 

radiation  more  attention  is  given  to  the  amplitude  function  A(k  ,k  ). 

x y 

The  MTF  for  an  optical  medium  is  simply  the  magnitude  function 

A(k  ,k  ) normalized  to  give  a value  of  unity  at  k = k = 0,  where  the 
x y x y 

object  signal  is  a spatial  sine  wave.  Definitions  for  MTF  and  related 
functions  are  given  below. 

Modulation  Transfer  Function  (MTF)  - a spatial  frequency  function 
which  gives  the  normalized  sine  wave  response  of  an  optical  element. 


Point  Spread  Function  (PSF)  - a spatial  function  which  describes 


how  an  optical  element  spreads  or  distorts  an  ideal  point  source.  May 
be  defined  as  the  normalized  intensity  distribution  in  the  image  plane 
of  an  ideal  point  source  in  the  object  plane. 

Edge  Function  (EF)  - the  normalized  response  of  an  optical  element 
to  an  optical  step  function  (edge).  EF  is  a spatial  function. 

Contrast  Transfer  Function  (CTF)  - a function  of  spatial  frequency 
which  gives  the  normalized  response  of  an  optical  element  to  a square- 
wave  input. 

Properties  of  the  Above  Functions 

a.  0 < MTF  < 1 

b.  0 < CTF  < 1 

c.  MTF  of  system  elements  can  be  cascaded  to  give  system  MTF. 

d.  CTF  cannot  be  cascaded 

e.  CTF  easier  to  measure  since  optical  square  wave  is  easier  to 
generate  than  optical  sine  wave. 

f.  CTF  can  be  converted  to  MTF. 

g.  EF  is  the  integral  of  PSF. 

h.  MTF  and  PSF  form  a Fourier  transform  pair. 

1.  EF  easier  to  measure  than  PSF  since  ideal  point  source  is 
impossible  to  model. 


APPENDIX  D:  TRANSMISSION  FUNCTION  FOR  A THIN  LENS 


If  t(x,y)  is  the  lens  transmission  function  assume 


t(x,y)  = p(x,y)  e 


j$(x,y) 


where:  p(x,y)  is  the  pupil  function  and  represents  the  aperture 

without  the  lens. 

+(x,y)  is  the  phase  function  which  depends  on  the  optical 


thickness  of  the  lens. 


SPw&f/c/U 

SujfMCfS 


R ' * xa4y*  t (*-*)' 
K*  * y*  + (RWO* 


OPTICAL 

AXIS 


- dMAV  - Tp 


Then  the  total  phase  del.,  ough  the  lens  is 


♦ * 2xknd(x,y)  *=  2*kd(x,y)  + 4(x,y) 


Now  the  additional  phase  delay  through  the  lens  Is  (over  what  would 
occur  without  lens) 

♦(x,y)  « 2wk(n-l)d(x,y) 


where:  d(x,y)  Is  given  in  (D2)  and,  as  can  be  checked  in  optics 

references,  we  have 


optical  path  length  * (refractive  index  • actual  length) 


Approximation:  for  a "thin"  lens 


Equation  (D6)  gives  thickness  for  parabolodial  surfaces  so  that  is  the 
approximation.  Using  (D6)  in  (D4)  gives 


If  define  focal  length 


Then  the  lens  transmission  function  is 


THIN  LENS  TRANSMISSION  FCN 


where  p(x,y)  is  the  pupil  function 


CONVENTION  ON  LENS  RADII  AND  FOCAL  LENGTH 


(DIO) 


where  both  R-j  and  R^  were  convex  surfaces  to  left  and  considered  positive 


for  surfaces  concave  to  left,  the  radii  should  be  negative 


ti£G*P*IE 


■ PfejCTM te 

The  conditions  under  which  the  focal  length  is  negative  can  be  found  by 
examining  eqn. (DIO)  eg. 


ffi 

*£«?/»  T)l>£ 


Foim'se 


Equation  (D9)  is  the  desired  result. 
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